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Introduction 



Computation is a physical process and the notion of a computable function, say /, rehes 
on the possibihty of implementing aphysical process transforming the input state in the 
desired output state /{Sq) [3, HJj. n_ 
The Church- Turing hypothesis [H) ^jB]- 

Every function that would naturally be regarded as computable can be computed 
by a universal Turing machine 

makes itself an implicit physical assertion which is explicitly stated in the Church- Turing- 
Deutsch hypothesis (a): 

Every finitely realizable physical system can be perfectly simulated by a universal 
model computing machine operating by finite means. 

Following this principle, a universal computing machine operating by finite means is able to 
simulate the evolution of a bundle of, say, n interacting electrons. Feynman showed \^\ that 
the complexity of the simulation, on a classical computer, of a quantum mechanical system 
scales exponentially with its dimension; but Feynman himself pointed out that quantum sys- 
tems are more "suited" for the simulation of other quantum systems; by "suited" he meant 
that a logarithmic reduction of the complexity of the simulation is achievable by means of a 
universal model quantum computing machine. 

Quantum computation origins from a question which is strongly suggested by Feynman's con- 
siderations: are there other hard computational problems which can exploit the features of 
quantum mechanical systems to be efficiently solved? 

Since the seminal work of Feynman quantum computation has known an enormous growth 
and nowadays it is a mature and vast research field in between physics and computer science. 
We refer to [lO| for an exhaustive introduction to the field of quantum computation, quantum 
information and quantum communication. 

Another problem motivating the investigation of quantum computing devices is the techno- 
logical advancement of semiconductor industry. A state of the art MOSFET (Metal Oxided 
Semiconductor Field Effect Transistors) has dimension of order lO^^m, the Bohr radius is 



^In f^l Toffoli illustrated the equivalence between a physical experiment and a computing process by means 
of the following suggestive example: 

Let us suppose that intelligent beings are observing us from a far away star. How could they 
understand whether we are carrying on a computation or a physical experiment? They could not 
understand it from what we are doing since there is no objective difference. The difference is in 
our intentions, in our knowledge, in our expectations. 
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approximately 10~^^m; following Moore's lawH, which has been, up to now, very accurately 
verified, in a few years the dimension of a MOSFET will reduce to fractions of a nanometer. 
At this scale, the evolution of the gates will be described by quantum mechanics. Indeed, 
quantum corrected diffusion models have already been introduced in the design process of 



nanoscale semiconductor devices (see Ul for an updated list of references) 



In this work we present the quantum mechanical computer proposed by Feynman in 1985 
and, since then, widely cited but seldom used. The main feature of the model is the presence 
of a builtin clocking mechanism managing for the ordered application of the computational 
primitives to the input/output register. 

In fact, given a transformation A to be applied to the input/output register, quantum com- 
putation starts from the decomposition of A into the sequential application of simpler com- 
putational primitives, Ui, U2, ■ ■ ■ , Un, which are unitary operators acting on few qubits of the 
input/output register at a time. The evolution of the input state into the output state is 
then seen as a discrete, stepwise, process: at any time step an operation is performed on the 
input/output register. 

In general, in quantum mechanics, the outgoing state at time t for a system with time inde- 
pendent Hamiltonian H is e~*^*| ipo ) , where | V'o ) is the input state; in other words, the state 
of the system evolves under the action of the unitary group generated by the Hamiltonian 
H. It appears to be very difficult to find, for a given special time i, the Hamiltonian which 
will produce A = e~*^* when ^ is a product of noncommuting matrices Ui,U2, ■ ■ ■ ,Un, from 
some simple property of the matrices themselves H. 

If one accepts time dependent Hamiltonians, it is fairly straightforward to write a Hamilto- 
nian for the evolution of the input state into the output state via the indicated intermediate 
steps; keeping in mind the mechanism of a synchronous system, one can imagine an inter- 
nal clock turning the interactions on and off. This idealized clocking mechanism is however 
not satisfactory: a classical macroscopic clock would destroy the coherence of the quantum 
system; so, to be consistent, the clock itself should be quantized; but as soon as we do it, it 



becomes clear that this clock would be affected by whatever it interacts with 13|] and time 
steps would be blurred. 

The problem of explicitly defining a time independent Hamiltonian driving the input state 
to the output state through the intermediate states determined b y th e ordered application of 
the computational primitives was overcome by Feynman in 1985 [14| . 

There are other interesting aspects of Feynman's proposal for a quantum computer. For 



example, it has been observed by Margolus 15|] that Feynman, in his model, 'managed to 



arrange for all the quantum uncertainty[...] to be concentrated in the time taken for the 
computation to be completed, rather than in the correctness of the answer'. 
More recently, Levitin and Margolus [3] related the maximum rate of information processing 
by a quantum computer to the available, conserved, energy. It is therefore of some theoretical 
interest to revisit a model, such as Feynman's, based on a closed system, evolving according 



^Moore's Law is the empirical observation that the transistor density of integrated circuits, with respect 
to minimum component cost, doubles every 24 months. It is attributed to Gordon E. Moore, a co-founder of 
Intel. 

^For example, following the method proposed by Benioff 12] it is possible to define a time independent 
Hamiltonian guiding the desired evolution from any input to any output state. The problem is that the explicit 
construction of such a Hamiltonian requires the prior knowledge of every step in the solution of every problem 
which the computer can solve 
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to a time independent and, therefore, conserved Hamiltonian 



Moreover, the doubt has been raised by Ahcki 17| that 'the idea that the physical time[...] of 



computation is proportional to the complexity [...] is only true of the existing digital comput- 
ers which are ensembles of controlled bistable elements which [...] can literally mimic logical 
operations'. Feynman's model provides an ideal context for the study of this issue: timing is 
modeled by a cursor, which jumps along a sequence of sites, indicating that the corresponding 
discrete operations should be applied. 

The thesis is organized as follows. 

In Chapter [T] we present the model, the basic clocking mechanism and establish our notation. 
In Chapter [5] we use Grover's algorithm study to introduce the full model, which we 

call an interacting XY system; particular attention is paid to the role of additional controlling 
spins in implementing successive visits to selected parts of the flow chart, or graph, of the 
algorithm in iterated computations (quantum subroutines) and to the locality of implementa- 
tion. 

In Chapter [3] we study the dynamics of Feynman's quantum computer, the timing and syn- 
chronization problems related to a rescaling of the clock to the quantum regime and propose 
a measurement scheme for the storage of results of computation via telomeric chains. 
In Chapter Owe pay specific attention to non-positional observables of the system: our main 
concerns will be speed (of computation), entropy (of controlled and/or controlling subsystem) 
and energy (of the system). In particular we relate the speed of computation to the group 
velocity of the cursor wave packet along the graph and discuss the buildup of entropy in 
the clocked subsystem caused by the spreading of the wave packet of the clocking agent. An 
outline of possible choices of the initial form of the wave packet bringing the entropy buildup 
close to a minimum is given. 

In Chapter [5] we consider the observable number of particles (agents performing a quantum 
walk along an XY spin chain), discuss the interest and limitations of the proposal of a multi- 
hand quantum clock (or multi-agent spin networks) as a substitute for the loops implementing 
iterated applications of quantum subroutines. 

The Conclusions and Outlook chapter is devoted to an exposition of open problems and future 
line of research. 
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Chapter 1 

The Feynman machine 



We present Feynman's model of a quantum computer stressing the role of the clocking mech- 
anism. We introduce the notion of logical successor of a given state and the related notion of 
Peres' constants of motion. We furthermore set the notation used throughout this work. 



1.1 The basic model 



It has been shown by Feynman [14| that it is possible to implement the sequential application, 
in the desired order, of the sequence 

Us-i- ...■U2-Ui = A 

of unitary operators to an input/output register by using s additional degrees of freedom: the 
program counter sites. 

For the sake of definiteness, we will think of each program counter site j = 1,2, ... ,s as 
occupied by a spin-1/2 system r(j) = (Ti(j), T2(j), r3(j)). We will refer to the collection of 
such spins, which act in effect as a quantum clocking mechanism, as to a program line. 
The input/output register will be, similarly, implemented by a collection of a certain number 
H of spin-1/2 systems o;(z) = ((Ti(i), o"2(i), cT3(i)), i = 1, 2, . . . , /i. 

We remind that the angular momentum operators r(-) satisfy the commutation rules of the 
Lie algebra on SU2 

[Th{x),Tj{y)] = 5^^yehjk^Tk{x) (1.1) 

ehjk being the Levi-Civita symbol and 5x,y the Kronecker delta. Obviously the same conditions 
are satisfied by the a's. 

We reflect the functional separation of the subsystems by calling Tlcursor 

= w®^with?^ = c^ 

the Hilbert space in which the collection of the spins making up the cursor are defined; 
analogously we refer to the Hilbert space the register is defined in as TCregister = 'H^'^. The 
overall system, register + program line evolves in the space TLmachine = '^register ® 'Hcursor 
under the action of the Hamiltonian 

^ = - 2 (E ® + I)^-(J') + ^i' ® + 1)) • (1-2) 

i=i 

where 

^^,^)^na)±Ma (1.3) 
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are respectively the raising and lowering (or excitation creation and annihilation) operators 
acting on the j-th. spin of the cursor and A is a scalar coupling constant (for notational 
convenience, we will set A = 1 unless otherwise specified). For the sake of definiteness, for 
every spin of the system we will take the eigenstates of the ts (respectively (73) operator as 
the basis for the Hilbcrt space H of a single spin and indicate them by the eigenvalue ±1 of rs 
they belong to. This basis is conventionally referred to as computational basis and the Pauli 
operators ti, T2, rs have the usual matrix representation: 

The evolution of the computing device is given by the solution of the Cauchy problem: 

Ui\m) = H\m) .15. 

llV'(o)) = |^i) ^'^ 

where we have set ^ = 1 for notational convenience. 
We define the operator number of excitation 

-^3 • '^machine ^ "^machine 

Ns = Ir^t'-^, (1.6) 
Ir being the identity on the register subspace. 

Proposition 1. The number of spins up, or excitations on the program line is a constant of 
motion; namely 

[H,N3]=0. (1.7) 

Proof. Since Ns acts only on the cursor, we can omit in this proof explicit reference to the 
register degrees of freedom. From the commutation rules defined above it follows: 

s-1 



[H,N,] = [^r+(x + l)r_(x) + T+(x)r_(x + l),^i±^] = 

x=l 

= 9 E E ([^+(^ + l)T_(x),r3(x)] + [T+(x)r_(x + l),r3(x)]) = 



2 

x=l y=l 
s-1 

2 



1 '"^ 

:;Y.i[r+{x + l)T^{x),T3{x)] + [T+ix+l)T.{x),T3{x+l)] + 



x=l 

+ [t+{x)T-{x + 1),T3{X)] + [t+{x)T-{x + 1),T3{X + 1)]) = 
1 

= 2^ir+{x + l)[r_(x), r^ix)] + r_(x)[r+(x + 1), r3(x + 1)] + 

x=l 

+ r_(x + l)[T+(x),r3(x)]+r+(x)[r_(x + l),r3(x+l)]) = 
= -r+(x + 1)t_(x) + r_(x)T+(x + 1) + 

+ r_(x + l)r+(a;)-r+(x)r_(x + l) = (1.8) 

□ 
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A first consequence of proposition [T] is that if the initial state of the cursor belongs to the 
N3 = k subspace of TLcursor , the state vector of the cursor remains, for every time t, in the 

subspace of Hcursor spanned by the ^ ^ ^ basis vector 

T+ijl)T+{j2) ■ . . T+{jk)\ ), 1 < jl < J2 < • • • < jfc < S (1.9) 

where with | ) we indicate the "all-down" state, that is the state 

|r3(l) = -l;T3(2) = -l;...;T3(s) = -l). 

The dimensionality of the state vector, and thus the complexity of any simulation of the 
system, is therefore significantly reduced. In particular, we will usually restrict our consider- 
ations to initial conditions belonging to the A''3 = 1 subspace and refer the cursor subspace 
to the orthonormal basis 

{|C(i))=T+(j)|0), l<j<s}. (1.10) 
The initial state of the computing device will be of the form 

|Vi) = |ii(i)))®|C(i)). (1.11) 

It helps the intuition to think of the initial state | i?(l) ) of the register as a simultaneous 
eigenstate of the components of the a spins in selected directions, encoding the initial word 
(or superposition of words) on which the machine is required to act. 

The intuition of "a single clocking excitation traveling along the program line" emerging 
from the above considerations is made precise by introducing the observable position of the 
excitation, or position of the cursor: 

Q = Ir^i23'-^^. (1.12) 

i=i 

Let us have a look on the Hamiltonian (jl.2p and on the interaction of the program line with 
the register degrees of freedom it describes. At any particular time t, if we expand e^**^ out 
as 

^-iHt h2^2^2 + ... (1.13) 

we find the operator H operating an innumerable arbitrary number of times and the total 
state of the system is a superposition of this possibilities. 

To illustrate the functioning of the clocking mechanism, we will consider only states of the 
form 

= \R{k))'^\C{k)), l<k<s, (1.14) 

where by | R{j) ) we indicate the {j — l)-th logical successor of the initial state | i?(l) ) of the 
register, that is 

I R{j) ) = . . . U2Ui\ R{1) ), (1.15) 

and the state of the cursor is a basis vector of the A''3 = 1 subspace of TCcursor ■ 

We begin with a look at the action of -ff on a state of the form (I1.14|) . with 1 < A; < s: 

s-l 

H\il^k) = ^C/x®T+(x + l)r_(x)| Vfc) + f/x"^®r+(x)T„(x + l)| V-fe) 
x=l 

= Uk\ R{k) )0\C{k + l)) + R{k) )0\C{k-l)) = 

= \R{k + l))(^\C{k + l)) + \R{k-l))(g)\C{k-l)) = 

= I Vfc+i ) + I ^fc-i ) (1.16) 
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We intentionally left out the case in which H acts on the states | ipi ) and \ ips)- In those 
cases, due to the boundary of the system, we have 

H\^Pi) = IV2) (1.17) 
H\i;s) = \i^s-i). (1.18) 

The coupling of the register with the cursor degrees of freedom seems to be of the following 
kind: if the excitation moves one step further, the logical state of the computation advances 
by one; if the excitation moves one step backward, the logical state regresses to the previous 
one. This property holds also when applying higher powers of the Hamiltonian operator to 
the system. For example, if we consider H^ and expand it we get 

s-l s-l 

= J]^^^j/^^+(^ + l)T-(:r)T+(y + l)r_(y) + 
x=l y=l 

+ U^U-^ r+(x + l)T_(x)T+(y)r_(y + 1) 
+ U~^Uy r+(x)r_(x + l)r+(y + l)T_(y) 

+ C/-i[/-i®r+(x)T_(x + l)r+(y)r_(y+l). (1.19) 

It is straightforward to see that, due to the commutation rules (jl.ip and to the conservation 
law ()1.7p . only some of the terms with x = y and x = y ± 1 survive. Thus, it is possible to 
simplify (|1.19|) getting 

s-l 

H^ = ^C/:,.C/^_i ®T+(x + l)r_(x)r+(x)T_(x- 1) + 

x=2 

= Ir T+{x)t_{x + 1)t+{x + 1)t^{x) + 

= r+(x + l)T_(x)T+(a;)T(X + 1) + 

= U-},U-},0t+{x-1)t^{x)t+{x)t,{x + 1). (1.20) 

For example, if H acts twice on the state | C{k) ) \ R{k) ), with 2<A;<s — 2we get 

H^\R{k))(^\Cik)) = Uk+2Uk+i\Rik))®\Cik + 2)) + 

+ 2\R{k))0\Cik)) + 

+ U,},U,'\R{k))0\C{k-2)) = 

= |V'fc+2) + 2| Vfc) + lV'fc-2). (1.21) 

Once more: if the position of the excitation is shifted by j-positions, j £ {0, 2} the logical 
state evolves or regresses accordingly by j-steps. This property extends to every power of H; 
in fact it has been shown by Peres [18] that, once defined the projection operator on the k-th 
logical state | R{k) ) 

Pk\R(,l)) = 6k,i\R{k)). (1.22) 

that satisfies 

Pk = UkPk-iU^^ (1.23) 

and the operator 

s 

P = ^Pk^\C{k)){C{k)\ (1.24) 

k=l 

the following holds 



1.1. THE BASIC MODEL 



5 



Theorem 1.1.1. 

[P,H]=0 (1.25) 

Before showing the proof, we observe that, given p.l4|) and proposition[Tl the Hamiltonian 
(II. 2p is equivalent to 



H = ^Ux(g)\x + l){x\ + U~^ (g)\x){x + l\ = 

x=l 
s-1 

= ^U^\x + l){x\ + U-^\x){x + l\ (1.26) 



x=l 



where | x ) is the eigenstate of the position operator Q, defined as in (I1.12p . belonging to the 
eigenvalue x, and the tensor product symbol has been dropped to shorten the expressions. 
Equivalently ()1.24|) can be rewritten as 

s 

P = Y,Pk\k){k\. (1.27) 

k=l 

Proof. We compute explicitly P H and H P. 

s s—1 

PH = Y,Pk\k){k\{Y,U^\x + l){x\ + U-^\x){x + l\) = 

k=l x=l 

s s—1 

= ^^PkUx\ k){k\x + l){x\+ PkU-^\ k){k\x){x + l\ 

k=l x=l 

s-1 

= J^P,+i[/,|x + l)(x|+P,C/-i|x)(x + l| (1.28) 



x=l 



s-1 



HP = {J2u^\x + l){x\ + U~^\x){x + l\)Y,Pk\k){k\ = 

x=l k=l 

s—1 s 

= ^^t/xPfe|x + l)(x| A:)(A;| +C/-i|x)(2; + l| A:)(A:| 
x=l k=l 
s-1 

= J^C/x.Px|^ + l)(x|+[/-ip,._i|x)(x + l| (1.29) 



x=l 



Thus 



PH-HP = {P,+iU,-U,P,)\x + l){x\ + iP,U-^-U~'P,+i)\x){x + l\ = 
= {U,P,U-'U, - U,P,)\ x + l){x\ + 

+ (P^U-' - U-'U,P,U-')\ x){x + l\=0 (1.30) 



□ 
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The space spanned by the initial state and its logical successors is, therefore, a constant of 
motion. The set {| ^fc ), 1 < < s} forms a complete orthogonal basis for the s-dimensional 
subspace of 'Hmachine effectively visited during the computation (the orthogonality of different 
basis vectors following immediately form ( C{j) \ C{k) ) = Sj^k)- It is worth mentioning here 
that the set of logical successors of the initial state, or Peres ' basis , can be algorithmically 
constructed; in fact, if we split the Hamiltonian (I1.26P into 



s-1 



H 



forward 



H 



backward 



Y,Ux\x + l){x\ 

x=l 

s-1 



x=l 



the set of logical successors of the initial state | ipi ) corresponds then to the set 

{lV'.) = ^)-LjV'i), l<A:<s} 
The Hamiltonian (|1.26|) can be rewritten using the Peres basis as 

s-1 



H = '^ \ tpk+i ){ipk\ + \ipk){ tpk+i 



(1.31) 
(1.32) 

(1.33) 

(1.34) 



k=l 



or as a s X s bi-diagonal matrix 



H 






1 


. 


. 





1 





1 . 


. 








1 


. 


. 











. 


. 


1 








. 


. 1 






(1.35) 



which is, up to constant diagonal terms, the finite difference approximation of the Laplace 
operator. The evolution of the system is thus of the form 



m)) = Yl k;s)\iJk) = Yl k; s) \ R{k) ) ® | Cik) 



(1.36) 



k=l 



k=l 



c{t, k; s) being a numerical functions of the time t, of the label of the logical successor k and 
parametric with respect to the length of the program line. For the sake of simplicity we 
postpone the discussion of the c{t, k; s) coefficients to chapter [3l 

From (|1.36|) it becomes clear how the clocking mechanism works: if a measurement is per- 
formed and the cursor is found at position N, then the register collapses into in the — 1 
logical successor of the initial state | R{1)). In particular, if the cursor is found in the last site, 
the s-th in our notation, the logical state of the register corresponds to the desired output 
state, that is A\ R{1) ). In Feynman's words (adapted to our notations), (jl.36p says that, 
starting from the initial condition (II. lip . "If at some later time the final site s is found to 
be in the |r3(s) = +1) state (and therefore all the others in |t3(j) = —1) ), then the register 
state has been multiplied by Us-i U2 • Ui as desired". 
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1.2 Continuous time quantum walks 

In this section we give an introductory overview on continuous time quantum walks (as 
op pos ed to discrete time, or coined, quantum walks which we will not deal with here). We refer 
to 19|,|20(] for an exhaustive treatment of quantum walks and of their algorithmic applications. 

Markov chains or random walks on graphs have proved to be a fundamental tool, with broad 
applications in various fields of mathematics, computer science and the natural sciences, 
such as mathematical modeling of physical systems, simulated annealing, and the Markov 
Chain Monte Carlo method. In the physical sciences they provide a fundamental model 
for the emergence of global properties from local interactions. In the algorithmic context, 
they provide a general paradigm for sampling and exploring an exponentially large set of 
combinatorial structures (such as matchings in a graph), by using a sequence of simple, local 
transitions. It is thus natural to ask whether quantum walks might be useful for quantum 



computation. In [2]|, for example, it is shown that there are graphs for which the time for 
a quantum walker to propagate between a particular pair of nodes is exponentially shorter 
than the analogous propagation time needed by a classical walker. 

A continuous time classical random walk on a graph is a Markov process. A graph is an 
ordered couple {V,E), where V is the set vertices, say {1,2,... and E a set of edges 
between vertices. 

A step in a classical random walk on a graph only occurs between two vertices connected by 
an edge. Let 7 denote the jumping rate. Starting at any vertex, the probability of jumping 
to any connected vertex in a time e is 7e (in the limit e — > 0). This random walk can be 
described by the s x s infinitesimal generator matrix M defined by 

Mab = < 



—7, if a 7^ 6, a and b connected by an edge 

0, if a 7^ 6, a and b not connected (1-37) 
^k^, if a = 6, A; is the valence of vertex a. 



lipa{t) denotes the probability of being at vertex a at time t, then it evolves under the master 
equation 

^ = -Y.M'^bMt). (1.38) 

b 

Following [22], a natural quantum analogue to the classical random walk described above is 
given by the quantum Hamiltonian with matrix elements 

{a\H\b) = Mab (1.39) 

I a) and | b) belonging to an assigned basis 1 1 ), | 2 ), . . . , | s ) of a f-dimensional Hilbert space. 
The Schrodinger equation for | tp^t) ) can be written as 

b 

We observe that, in some sense, any evolution in a finite-dimensional Hilbert space can be 

thought of as an oriented graph with Hermitian weights. 

We furthermore point out that whereas (11.380 conserves the probability 
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the Schrodinger equation (ll.40p preserves probability as the sum of the amphtudes squared 

Y.\{a\m)\'- (1-42) 



The simplest graph we can take into account is the one-dimensional lattice Z, resulting in a 
nearest neighbor Hamiltonian defined by 



2| J ) + 1 i - 1 )), 



(1.43) 



that is the discrete approximation of the Laplace operator dP/dx'^, with A = y^2/A, A = 27. 
In the basis |l),|2),...,|u), the matrix representation of the Hamiltonian specified in p.43p 
is 

/ -2 1 ... \ 
1 -2 1 ... 
1 -2 ... 



H 



1 





V 













1 



(1.44) 



-2 / 



which is, up to an additive constant, equivalent to the Hamiltonian p.35|) of the basic model 
of the Feynman machine. We point out that in ()1.35p the basis states were the states \ipk), 
with 1 < k < s, defined as in (|1.33p . namely the Peres basis. 

We delay to chapter 3 the discussion of the dynamics of an excitation traveling through the 
linear chain or, more generally, on graphs of the form of figure [231 



Summary 

We have presented the clocking mechanism of Feynman's quantum computer: the numerical 
coupling constants between nearest neighbor sites of a spin chain with XY interactions are 
substituted by unitary operators acting on the register subspace. The so created entanglement 
between the clocking and the register degrees of freedom implements the same kind of timing 
used in classical computers to manage the ordered application of computational primitives to 
the input/output register. 

We have made evident by means of the Peres basis, that the evolution of the clocking sub- 
system is independent of the operations performed on the register, as long as we use unitary 
operators on the latter. 

The continuous time quantum walks paradigm has been deeply investigated in recent years, 
with the hope of getting some new technique to define new quantum algorithm working faster 
than any classical one. The kind of interaction used in Feynman quantum computer extends 
the continuous time quantum walk: the walker does something while traversing the graph. 

In the following chapter we will extend the linear chain to more general graphs. This will 
allow the implementation of basic flow control mechanisms, such as the IF. ..THEN. ..ELSE 
and the iteration of quantum subroutines, and a simplification of the interactions between 
elements of the system. 



Chapter 2 

Simplifying the implementation 



As a case study, we present an implementation of Grover's algorithm in the framework of 
Feynman's cursor model of a quantum computer. Using the cursor degrees of freedom as a 
quantum clocking mechanism allows Grover's algorithm to be performed using a single, time 
independent Hamiltonian. We examine issues of locality and resource usage in implementing 
such a Hamiltonian. In the familiar language of Heisenberg spin-spin coupling on a linear 
chain of spins, we introduce occasional controlled jumps that allow for motion on a planar 
graph: in this sense our model implements the idea of "timing" a quantum algorithm using 
a continuous-time quantum walk. In this context we examine some consequences of the 
entanglement between the states of the input / output register and the states of the quantum 
clock. 



2.1 Grover's algorithm on a Feynman machine 

The star ting point of our discussion is the analysis of the physical aspects of Grover's algorithm 
given in [23t | and 22l |. 

Suppose one is given an "oracle" able to compute, in a quantum reversible way, the indicator 
function of a binary word a G {—1, l}'^ of an assigned length fi. We will assume, for the sake 
of definiteness, that this computation is performed by applying a unitary transformation A 
to the input /output register of length i/ = fi + l. Suppose that A results from the action, for 
a fixed amount t of time, of a Hamiltonian K{a.), that is: 

A = exp{-itK{a)). (2.1) 

It is then possible to arrange things in such a way that the state 

I a) = I 0-3(1) = ai, 0-3(2) = a2, . . . ,o-3(;u) = ) (2.2) 

that corresponds to having the word a written on the register, is the jjround state of K{a). 
The search for the ground state of K(a) is performed, in Reference [22], following the simple 
idea of perturbing the Hamiltonian K(aL), 

K{a) ^ K{a) + p (2.3) 

with a perturbation /3 chosen in such a way that a suitable initial condition oscillates about 
the state | a ) with a period proportional to 2^^"^ , becoming, at a time 0(2^/2)^ parallel to 
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the target state. 

By applying Trotter's product formula to (j2.3p . it is shown, in Reference 23|, that no sig- 
nificant loss in the probability P(a) of finding a, at suitable values of time t, results from 
alternating intervals of time in which only the "oracle" Hamiltonian K{a.) is active, thus in 
fact applying the "oracle" transformation A, with intervals in which only /? is active, thus in 
fact applying the "estimator" transformation 

B = exp(-it/3) (2.4) 

The oscillatory nature of the quantum search algorithm is confirmed, in this discrete time 
setting, by the analysis of Reference [24 1 . 

We have seen that the administration, in the correct order, of the oracle and estimation 
transformation A and B to the register can be realized by means of the clocking mechanism. 
In our implementation we define the operator A : TLregister — > "^register through the action 
on simultaneous eigenstates of o"3(l), cr3(2), . . . , a^[fi),a^{u): 



0-3(1) = 2:1,0-3(2) = Z2,.. . ,0-3(/i) = Zfj,,a^{u) = Zy) = 

J I 0-3(1) = ^1,0-3(2) = 22, • • • ,o"3(/i) = z^,aj,{v) = -Zy ), if z = a 
1 I 0-3(1) = 21,0-3(2) = 22, • • • ,o-3(/i) = 2^,0-3(1^) = Zu), if a 



(2.5) 



where z = (21, 22, . . . , 2^) G {-1,1}^, Zy £ {-1,1}. 

The oracle A performs a quantum reversible computation of the indicator function of a by 
flipping the component 3 of the output qubit o;(i/) iff the word a is written on the register in 
terms of the ct(1),oi(2), . . . ,a(/i) components of the input qubits. 

Define, in a similar way, a linear operator B : TCregister — *■ T~(-register through the following 
action on the simultaneous eigenstates of o-i(l), o-i(2), . . . , o-i(;u), 0-3(1/): 



B I 0-1(1) = Xl, 0-1(2) = X2, . . . ,0-l(/i) = X^, 0-3(1/) = Zy) = 

rrn { tA = — "V" = 1 

(2.6) 



o-i(l) = Xl, 0-1(2) = X2, . . . ,o-i(;u) = x^, 0-3(7/) = -Zy ), i/ X = 1^ 



[I 0-1(1) = Xl, 0-1(2) = X2, . . . ,o-i(/u) = x^, 0-3(1/) = Zy), if :x.^ 1^ 
where X = (xi, X2, . . . , x^) G {—1,1}'^, Zy £ {—1,1} and 1^ = (1,1,..., 1). Following the 

fi times 

prescription of section [TTT] we define the Hamiltonian 

s-l 

H = Y,Uj\j + l){j\ + h.c. (2.7) 

x=l 

where 

{A if j is odd 
■' ■' (2.8) 
B if j is even. 

We set the initial state of the machine to 

IVi) = \Rg{1))\C{1)) = 

= |ai(l) = l,ai(2) = l,...,ai(/i) = l,a3(z/) = -l)|Q = l), (2.9) 

where C(l) is defied by (jl.lOp and \ Q = 1) is the eigenstate of the position operator Q defined 
in p.l2p belonging to the eigenvalue 1. 



2. 1 . GROVER 'S ALGORITHM ON A FEYNMAN MACHINE 
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A 



s-1 



Figure 2.1: The Hamiltonian 12. 7[ ) describes an XY interaction between nearest neighbor cursor spins. 
A is the "coupling constant" between spins corresponding to odd links; B is the "coupling constant" 
between spins corresponding to even links. Both A and B are, in fact, functions of the register spins. 

Following (|1.33p we define the basis states \ ipj ), 1 < i < s of the subspace of 'Hmachine the 
system evolves in. We recall that the basis states are of the form 

\^k) = Uk-i-...- Ui\ Rg{1) )\Q = k)- (2.10) 

if k is an odd number, = 2n + 1, it is 

Uk-i ...Ui = {B- Af. (2.11) 

An explicit expression for {B ■ A)" | cjjfl) = xi, o"i(2) = X2, • • • , cri{fi) = ) can be found by 
the iterative procedure of Reference [2j] : 

(5-^)"|cti(1) =Xl,CJl(2) =X2,...,CTl(/u) =x^) = 

= (a„(M)|a)3 + /3„(/i) Ez^al^a)) (2.12) 

where 

a„(/i) = (-l)"sin((2n + l)x(^)) (2.13) 

= arcsin(2-''/2) (2.14) 

/3„(^) = -A=^cos((2n + l)x(/i)) (2.15) 

In (I2.12P we have omitted explicit reference to | cri{iy) = — 1 ), as the conservation law 

[H,aiii^)]=0. (2.16) 
allows us to do, and we have set, for every z = {z{l), z{2) . . . , z{fi)) £ {—1, 1}^ 

I z)3 = I 0-3(1) = zi,a3{2) = Z2,... ,cr3{n) = z^ ). (2.17) 
The case of an even value of /c = 2n + 2, can be similarly analyzed, by observing that 
A - (S- A)"|ai(l) =xi,cTi(2) = X2,...,ai(^) =x^) = 

= (-«„(/.)! a )3 + /3n(^) Ez^al^s)) (2-18) 
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Observation 1. The state of the register evolves in a proper two dimensional subspace of 
'^register ■ The evolution of Grover's register state can be simulated using a single qubit rotat- 
ing in a two dimensional Hilbert space 7i. We will exploit this feature later. 

Summarizing, and considering, for the sake of definiteness, the case of an odd value of 
s = 2g + 1 the sohition of the Schrodinger equation 

^J^\m)=H\m), (2.19) 
under the initial condition (|2.9|) can be written as: 

\^{-t)) = ^c(t,2n + l;s) ja„(^)|a)3+/?„(^) ^|Z3) I 

n=0 y z^a / 

® \ai{v) = -I) ®\Q = 2n + l) + 
( 

+ c{t, 2n + 2; s) -a„(^)| a )3 + ^ | Z3 

n=0 \ z^a 

(g, \ai{iy) = -1) (g>\Q = 2n + 2) (2.20) 

The conditional probability of reading upon measurement the word a on the register given 
that the cursor has been found at site k is then given by |Q:fc(/i)p. 



2.2 Iteration of quantum subroutines 

Grover's algorithm provides a quadratic speedup with respect to any classical algorithm for 
the search of a keyword in an unstructured database of 2^^ words [2^ : it is sufficient to apply 
the B ■ A operator 0(2^/^) times to the register to have the word a we are looking for written 
on it. Thus, if we sticked to the linear chain model, we would need an exponential number of 
functional blocks A and B and program line sites to perform the algorithm. 
In this section we show that this cost in terms of space can be made linear in fi by using 
quantum subroutines. 

For every non negative integer K, we wish to show a quantum clocking mechanism able to 
apply 2^ times the transformation BA to the register qubits g_{l),a{2) . . . ,a{v), by repeatedly 
using the same "piece of hardware" that applies BA just once. We will show that this clocking 
mechanism will involve 

s{K) = AK + 3 (2.21) 

cursor qubits T(l),r(2) . . . ,t{s{K)). 

In order to keep track of the progress of the 2^ executions of the assigned subroutine BA, 
there must be a subsystem (the subroutine counter) having 2^ different states: it will be 
constructed in terms of K qubits p{l),p{2) . . . ,p{K), p{j) = {pi{j), P2{j), PsU)), l<j<K- 
We will denote with TLcounter the 2^ dimensional state space of the counter degrees of freedom. 
The definition of the Hamiltonian operator on TLregister ® '^counter ® Ti-cursor will be given by 
an iterative scheme. Set, for i = 1,2, . . . , s{k) — 2 

hoih * + 2) = (^/oruiard(^) + B fc^^ard{i "I" 1)) ® Icounter 

= {A\i + l){i\+B\i + 2){i + l\)(^hounter (2.22) 




Figure 2.2: Do BA twice. 



where Icounter is the identity operator in Hcounter- The operator /io(i, i + 2) apphes the trans- 
formation BA to the register while the cursor jumps from site i to site i + 2 (see figure [T2]) . 
Fori = l,2,...,s(K)-6 

hi{i,i + 6) = p+(l)| i + I + 

+ p^{l)\i + 2){i + l\ + 
+ hoii + 2,i + 4) + 
+ />_(!)! i + 6)(i + 4| + 
+ p+(l)| z + 5)(i + 4| + 

+ p_(l)|i + l)(i + 5| (2.23) 

The term /9_(l)|i-|-6)(i-|-4|-|-/9+(l)|i-l-5)(i-|-4| in (j2.23p is an example of the implementation of 
a conditional jump through the SWITCH primitive. The first addendum acts non-vanishingly 
only in the subspace belonging to the eigenvalue +1 of the controlling qubit P3(l) and sends the 
excitation of the cursor from i-|-4 to i-|-6; the second addendum, in turn, acts non-vanishingly 
only in the subspace belonging to the eigenvalue —1 of P3{1) and sends the excitation of the 
cursor from z -|- 4 to i -|- 5. Notice that in this implementation of the IF. ..THEN, ELSE... 
construct, the controlling bit P3{1) gets inverted. 
The iteration step from /ij-i to hj is given by: 

hj{i,i + 4j + 2) = p+{j)\i + l){i\ + 

+ p^{j)\i + 2){i + l\ + 

+ hj{i + 2,i + 4j) 

+ p.{j)\i + 4j + 2){i + Aj\ + 

+ p+ij)\i + ij + + + 

+ p^{j)\i + l){i + Aj + l\ (2.24) 

and is represented in figure 12. 4[ For a fixed value of the positive integer K we define the 
forward part of the Hamiltonian as 

HforwardiK) = /^^(l, s(i^)) = /ij^ (1, 4E: + 3) (2.25) 
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Figure 2.3: Do BA twice. Here the operators acting on the counter degrees of freedom are encoded 
in the following graphical convention: to a raising operator corresponds a thick solid black line; to 
a lowering operator a dashed line; to the negation pi a thick solid green line. The functional block 
ho{i + 2, i + 4) is indicated as a coupling constant on the corresponding edge. 



and the Hamiltonian as 

H{K) = Hforward{K) + Hhackward{K) = Hforward{K) + Hforward{K)K (2.26) 

Let us consider the Schrodinger equation 

ij^\m) = -\H{K)\m) (2.27) 

and the initial condition 

|V(0)) = |ai(l) = l,cTi(2) = l,...,cTi(;u) = l,ai(z.) = -l)® 

® I />3(1) = -1, P3(2) = -1, . . . , P3 W = -1 > 

|<9 = 1). (2.28) 

Equation ()2.27p under initial condition ()2.28p is extremely easy to solve because of the con- 
servation laws 

[H{K),ai{v)] = (2.29) 

[H{K),N^] = (2.30) 

[H{K),P{K)] = (2.31) 

where the operator P{K) is the projector operator on the subspace of TLregister ® '^counter ® 
'Hcursor Spanned by the 2^+^ — 5 orthonormal vectors defined, similarly to ()1.33p . by 

{\^k) = H'f-i^J ^Pl),l<k< p{K)} (2.32) 

with I ^/^i ) = I 1^(0) ). Because of the above considerations, the solution of ()2.27p . p.28p will 
be of the form 

piK) 

l^(i)>= j;c(t,j;p(K))| Vi) (2.33) 
i=i 



2.2. ITERATION OF QUANTUM SUBROUTINES 
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Figure 2.4-' Do BA 2^ times while the cursor moves from i to i + s{j) — 1 

A full understanding of the solution ()2.32l) requires the analysis of the states | V'i )) for 
j = l,2,...,p{K). 

Because of (I2.29P all of them are eigenstates of cri{u) belonging to the eigenvalue —1; from 
now on, we omit the explicit reference to this fact, using the shorthand notation 

I o-i(l) = xi, . . .,ai{fj,) = x^,ai{u) = -1 ) = | o-i(l) = xi, . . .,ai{fi) = ) = | ) (2.34) 

for x^ G {-l,ir. 

Each of the vectors | 'i/'j ) will be, furthermore, a simultaneous eigenvector of each of the op- 
erators p3 = (/J3(l), p3(2) . . . ,p^{K)). Calling rj G { — 1, 1}''^ the collection of the eigenvalues 
to which I ipj ) belongs, we will write, for istance, 

|V'i) = |li)0|/O3 = -l)®|Q = l>- (2.35) 

and 

\ipj)= A'^ (5^)"^" \li)0\ps = rj)0\Q = qj). (2.36) 

The explicit iterative algorithm by which Cj, rij, qj, rj can be computed is strictly parallel 
to the iteration procedure of figure [231 For the following discussion it is sufficient to define 
the exponents ej and rij. 

I otherwise 



where, for i = 1,2, . 



2ir + 2 + 5(i-l) + 3^e2(x) = 

x=l 
K-l 

2E: + 2 + 5(i - 1) + 3 ^ ^(i - l)/2^-'' . (2.38) 



h=l 
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In (I2.38P we have indicated by 62(0;) the exponent of the prime factor 2 in the factorization 
of the positive integer x, and by [y\ the integer part of the positive real number y. 
Let us focus our attention on the states | \ tp . . . , \ ipj^,. ) ), 

IV'.i) = ^|li)|/>3 = li)|Q = 2K + 2) = 

= \U)\N,, = l)\Q = n). (2.39) 

In (j2.39p we have given a numerical meaning to the content of a subroutine counter by defining 
the operator 

Ar^^ = l + ^i±|3M2?/-i. (2.40) 

y=l 

In all the predecessors \ tpi),\ tp2) , ■ ■ ■ , \ V'ji-i )> I V'ii ) the register is in its initial state | 1 )i. 
The immediate successor of | V'ji ) is 

I V'.i+i ) = BA\l)i\Np, = l)\Q = 2K + 3). (2.41) 

In all of the states | i/jji+i ))•••)! "0^2-1 ) register remains in the state BA\ 1 )i ; the content 
of the register changes only at step j2 , where it is 

\tPj, ) = ABA\l)i\Np,^ =2)\Q = 2K + 7). (2.42) 

At each of the steps ji the state of the register gets acted upon by an additional A and at 
step ji + 1 by an additional B. In steps from ji + 2 to jj+i — 1 the state of the register remains 
unaltered. 

The content of the register becomes (BA)"^ for the first time at step + 1 = p{K) — K 
and such remains until the last step p{K). 

The exponent rij in ()2.36p is therefore equal to the number of "non-trivial" steps ji that 
precede step j: 

n, = |{l<i<2^:i, <i}|. (2.43) 

It is, therefore 

Uj = 0, for j <2K + 2 

n2K+3 = Ij 

Uj = 2^ , ior j > p{K) - K. (2.44) 
For 2i^+3 < j < p{K) — K, rij grows in an approximately linear way because of the inequality 

ji > 2K + 2 + 5{i-l) + 3(^l-^^{i-l)-3L{i) (2.45) 

< 2i^ + 2 + 5(z-l) + 3(^l-^) (i-1) (2.46) 

L{i) = Llog2(i-l)J (2.47) 

which easily follows from (I2.38h and from the fact that x — 1 < [x\ < x. 
This justifies the approximation 

'0, for 1 < j < 2if + 2 

l + p(4-3^-3 0'-(2^ + 3)). ioT2K + 3<j<p{K)-K (2.48) 
2^, for p{K) - K < j <p{K). 



2.3. EQUIVALENT 'LOCAL' HAMILTONIAN 
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The iteration of quantum subroutines mechanism provides thus a way to reduce the space 
complexity of an algorithm; this, in turn, reduces the complexity of the physical implemen- 
tation of the Grover algorithm on a Feynman computer. In the next section we address the 
problem of further simplify the implementation of quantum functional blocks. 



2.3 Equivalent 'local' Hamiltonian 

We restrict the class of unitary operators Uj acting on the register to rotation operators acting 
on a single spin. In this way each addendum in the Hamiltonian involves at most 3 bodies: 
two spins of the program counter and one of the register. The fulfillment of this requirement 
makes the architecture of the quantum computer modular [3] > thus simplifying the physical 
implementation of the desired interactions. We will show that this requirement does not affect 
the computational power of the Feynman machine. Furthermore, we will provide an analysis 
of the computational cost, in terms of space (additional qubits) involved in substituting such 
non-local terms with equivalent terms in which only interactions between two cursor spins 
and at most one register spin appear. 

The explicit expression, in terms of the register spins, of the oracle operator A defined in 
(|23]) is: 

A=l + iaA^)-l)fl'-±^. (2.49) 

j=i 

The analogous expression for the estimator operator B defined in ()2.6p is 

B = l + ia,i.)-l)fl'-±^. (2.50) 

i=l 

where 1 = is the identity operator in TLregister ■ 

In the Hamiltonian H defined in (12. 7p and in the Hamiltonian H{K) (I2.26p . there are, there- 
fore, non local terms such as ^| j -|- 1 )( j | and -B| j -|- 1 )( j | involving many-body interactions 
among two cursor spins and all the register spins. Our goal, in this section, is to show that 
both the transformation A and B can be implemented by means of local operations involving 
at most three bodies: two degrees of freedom of the clock and one of the register. 

For the sake of definiteness we concentrate our attention, to start with, on the clocked imple- 
mentation of the C — NOT logical operator 

CNOT{j, j + l) = C- NOT I J + 1 ) ( J I + ^-c- (2.51) 

where 

C-Aror= l + ((7i(i/)-l)f[ii|^. (2.52) 

i=l 

The C-NOT is a logical reversible binary operator whose truth table is given in table \2.1i\ . 
The action of CNOT{j, j + 1) on the computational basis is "Flip the component 3 of the ly-th 
qubit iff (73(1) input qubit points in the -|-1 direction, starting with the cursor in position j". 



The case /x = 1 of one controlling qubit has been studied in [IJ]. It involves the introduction 



^The C-NOT can also be seen as a reversible version of the binary sum. 
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/ -1 \ 

f^3(l) 




^3(1)' 


f^3(z^) 


-1 


-1 


-1 


-1 


-1 


+ 1 


-1 


+ 1 


+ 1 


-1 


+ 1 


+ 1 


+ 1 


+ 1 


+ 1 


-1 



Table 2.1: The truth table of the C-NOT; cr3(l) and 173(1') are the input values whereas 173(1)' and 
u^(v)' are the output values. 



C-Not 



o-(l) o-(l) 
j+1 (Ti(v) j+2 




Figure 2.5: C^NOT{j,j + 5); we have used the graphical conventions of fioure \2.3[ 



C-Nol 




Figure 2.6: C^'-^NOT CNOT. 



of si = 6 cursor qubits T_{j),T_{j + 1), . . . ,r(j + 5) and, supposing that the controlling qubit 
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is (73(1) and the controlled one is (73(1'), of the local Hamiltonian 

HcNOT 



— 






+ 


0-1 


H|j + 2)(j + l| + 


+ 




(l)li + 5)(j + 2| + 


+ 




(l)|i + 3)(j| + 


+ 


1 J 


+ 4)(i + 3| + 






delay line 


+ 




(l)li + 5)(j + 4| + 



z. (2.53) 

A graphical representation of ()2.53p is given in figure [231 

The term R(j + 3,i + 4) = 1\ j + 4){j + 3\ in (f233|) . represented as R(l) in figure [231 
plays the role of a delay line of length 1. It makes the length Ti = 4 of the computation 
independent of the input word in the sense that an initial state of the form | (T3(1) = +1 ) ® 
I •73(1^) = ) \ Q = j ) has the same number of logical successors 

\ a3{l) = -1) \ a-siu) = z^)\ Q = j + 1) 
I C73(l) = -1 ) » I ^3(1.) = -z,)\Q=j + 2) 

I c73(l) = +l)0\a3{iy) = -Zu)\Q=j + 5) (2.54) 

as an initial state of the form | (T3(1) = —1 ) | <73(j^) = Zu)\ Q = j ) , which has the successors 

I ^3(1) = +1 ) I asii^) = z^)\Q = j + 3) 
\ asil) = +1) \ a3{iy) = z^)\ Q = j + 4) 

I (73(1) = -1 ) I asiu) = z,)\Q=j + 5). (2.55) 

Figure [231 shows the iteration step leading from C^^^NOT to C^NOT through the intro- 
duction of the additional controlling qubit (J3 (/i) . The length of each computation increases 
from the previous value to 

r^ = T^_i + 2 = 2(^ + l). (2.56) 

The number of cursor qubits increases, because also of the delay line i?(T^_i), from the 
previous value s^-i to 

= 2 + s^_i + T^_i = (/i + 1)(^ + 2). (2.57) 
The iteration step C^'^NOT C^'NOT is exphcitly given by 

C'NOTijJ + 8^ - 1) = a_(/i)| J + 1 )( j I + C^^-'NOTij + + s^_i) + 

+ Cr+(M)I j + S^^ - + I + 

+ a+(//)| j + s^_i + l)(i| + 

+ Yl \j + s^-i + k + l){j + s^-i + k\ + 

k=l 

+ \j + s^-l){j + s^-l + T^^i\+H.c. (2.58) 
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Observation 2. For fi = 2 we implement the CCNOT gate, or Toffoli gate, which is a 
complete logical basis for the reversible boolean functions. The Feynman machine is thus able 
to compute at least every function computable by a deterministic Turing machine (or, in 
other words, the class of functions reversibly computable by a Feynman machine is at least as 
wide as the class of partial recursive functions ) by means of only three body interactions . In 
fact, every term in the Hamiltonian for the CCNOT will involve two spins of the program 
line and one spin of the register. The truth table of the CCNOT is reported in table \2.^ 



^3(1) 


^3(2) 




^3(iy 


^3(2)' 
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Table 2.2: The truth table of the C-C-NOT. 




The "double diamond" circuit of figure 12.71 is an example of a planar graph implement- 
ing, with the right setting of the input/output register, Grover's algorithm for fj, = 3: the 
excitation, which in our setup will initially be located at site 1, moves on the planar graph 
depending on the state of the input /output register. During the "walk", the excitation inter- 
acts with the register degrees of freedom. The implementation of an algorithm on a Feynman 
quantum computer thus leads to the definition of interacting quantum walks. 
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2.4 The SWITCH and the projective control 

Although the CCNOT is universal with respect to the class of reversible boolean functions, 
it is useful to generalize the kind of trajectory control used for the CCNOT to the "pure" 
IF. ..THEN. ..ELSE construct. Consider the Hamiltonian 

HswiTCH = CT_| j + l)(j|+A| j + 2)(j + l|+a+| j + 5)(j + 2| + 

+ fT+|j + 3)(j|+5|i + 4)(i + 3|+(j+|i + 5)(j + 4|+/i.c. (2.59) 

It is worth mentioning that the same control on the trajectory of the walker can be imple- 



Switch 




j+2 



\ 



j+4 



Figure 2.8: A SWITCH implementing the IF{a = +1) THEN A ELSE B instruction of Hamiltonian 
\2.59\) . A and B are primitive acting on the register. It is required that A and B do not modify the 
state of the controlling qubit a. 

mented by means of a purely projective mechanism; consider the following Hamiltonian 



HsWITCH 



+ 



+ 





2 


1 


+ 0-3 




2 


1 


- 0-3 




2 


1 


- 0-3 



|i + l)(j|+^| j + 2)(j + l| + 
i + 5)(i + 2| + 
j + 3)(j|+S|j + 4)(j + 3| + 
i + 5)( j + 4| +h.c. 



(2.60) 



Although the control of the trajectory of the cursor traveling along the program line is the 
same as the one of (I2.59p . the state of the controlling qubit a is left unchanged all through 
the SWITCH. This projective control mechanism has been introduced in |25| and makes it 
easier to track different computational paths; for example, to understand what happens on the 
upper branch of a switch it suffices to consider only the projection on the CT3 = +1 subspace 
of the state space of the machine. 
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Summary 

Throughout this section we have used a top-down approach to the implementation of Grover's 
algorithm on a Feynman machine: we started from the computational primitives A and B seen 
as black boxes and ended up with the same primitives decomposed in three body interactions. 
Moreover, we have shown that the iteration of quantum subroutines can be implemented by 
means of a linear amount of space resources. The possibility of using the same functional 
block many times amply justifies this linear space, an therefore time, overhead. 
The topology of a Feynman machine can then be extended from a linear graph to a planar 
graph; at the same time, the interaction between the register and the cursor can be exploited 
to make the clocking excitation visit different computational paths; the possibility of having 
the controlling qubit in a superposition of states allows for the different computational paths to 
be explored simultaneously. In this framework, the meaning of quantum parallelism becomes 
clear: different physical computational trajectories are visited simultaneously. 
The algorithmic construction of the Peres basis is directly extensible to this general form of 
graph. This allows for the state of the system to be described by (|2.33p . 



Chapter 3 



Quantum timing and 
synchronization problems 



The interacting quantum walk presented so far implements a clocking mechanism. In fact "If 
at some later time the final site s is found to be in the | C(s) ) state then all the computational 



primitives Ug-i U2-Ui have been administered to the register in the correct order" [IJ]. 

We have shown in chapter [2] that the same kind of interaction can be exploited to control 
the motion of the clocking excitation on a planar graph. However, the quantum nature of 
the cursor imposes limitations on our ability to know, without preforming a measurement, 
whether the computation has finished. In this chapter we study the motion of the cursor and 
make quantitative the following qualitative assertions: 



(i) at no instant of time the probability \c{t, s; s)p is larger than const ■ s 



(ii) the cursor keeps bouncing back and forth between positions 1 and s, thus in effect 
making the above upper bound attainable only at selected instants of time. 

In other words, by scaling down the clocking mechanism of the computing device to the 
quantum regime two quantum phenomena become relevant: the spreading of the probability 
distribution of the excitation (or pointer) along the program lines, and the scattering of the 
probability amplitude at the two endpoints of the physical space allowed for its motion. 



We begin this chapter by discussing the timing problems (i) and (ii) We then propose a 



measurement scheme which, as proved in |26l ]. makes the upper bound on the probability cost 
of the implementation of an algorithm on a Feynman machine less severe. 
We conclude the chapter by showing that the results obtained for the linear chain can be 
extended to a more general class of planar graphs, such as the one of section [2l by means of 
synchronizing delay lines. 

3.1 Motion of the cursor 

The eigenvalue problem for the Hamiltonian 

A 

H = --Y,Ucc\x+l){x\+h.c., (3.1) 

x=l 
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Figure 3.1: We take a linear chain of length s — 100. Figures (a) (h) and (c) represent the probability 
distribution for the position observable Q at different times; the fast spreading of the wave packet is 



confirmed by the; 3.1(a) 3.1(b) 3.1(c) represent the probability of finding the clocking excitation at 



sites 1, 50 and 100 respectively. The reader should pay attention to the vertical axes scale. 



corresponding to a sequential program line of length s, is solved by the following Ansatz for 
the eigenstates: 

s 

\e) = Y,v{x)U^-iU^...U2Ui\r)(^\C{x)), (3.2) 

x=l 

where | r ) is a unit vector in TCregister and | C{x) ) is defined as in (ll.lOj) . together with the 
boundary conditions f (0) = v{s + 1) = 0. Inserting this Ansatz, the eigenvalue problem 

H\e) = e\e) (3.3) 

becomes: 

e v{x) = --{v{x + 1) +v{x - 1)). (3.4) 



This leads in an obvious way [27[ to the eigenvalues 

Cfc = -A cos ^^^^^ , A: = 1,2,... ,s. (3.5) 

The multiplicity of each eigenvalue is equal to d = dim{7iregister )• An orthonormal basis in 
the eigenspace belonging to the eigenvalue e/^ is given by: 

s 

I '^k;rj ) = ^Vk{x)Ux-i • . . . • ) (g) I C{x) ), (3.6) 

x=l 

where \ ri) , . . . , \ r^) is an orthonormal basis in Tiregister , and 
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We observe that, for 1 < xq < s, it holds 

d 



C(xo) 



s + l 



■ sm 



/ k7rxo\ 
\s + l) 



For xq = 1 we obtain for the coefficients c{t, k; s) 



c{t, x; s) 



s + 



k=l 



iXt cos 



kn 
s + l 



sm 



kn 
s + l 



(3.8) 



sm 



kirx 
s + l 



(3.9) 



We point out that, given an Hamiltonian of the form ()3.ip and an initial state \ ipi ) = 
I R{1) ) ® I C{1) ), once defined the corresponding Peres' basis {tpk}k=iy the dynamics of the 
system is always of the form 



^(*) ) = Yc{t,x;s)\ V'x ), 



(3.10) 



x=l 



thus independent of the transformation applied to the register. 

This property of the computational device allows for the timing problem to be discussed 
independently of the transformations carried on the register, as long as we act with unitary 
operators on the latter. 

With the analysis of Reference ^] it has been shown that: 

Theorem 3.1.1. Given a linear chain of length s evolving under the Hamiltonian H of il.35\) . 

the probability P{t,s;s) of finding, upon measurement, the clocking excitation at the final site 
s is bounded by 

P{t, s; s) = \c{t, s; s)\'^ = 0(s~t). (3.11) 
We refer to Reference ^2^ 1 for a proof of the statement. 

Theorem 13.1.11 makes it clear that, as soon as the clocking mechanism of a quantum 
computer is scaled down to the quantum regime, a new problem, to which we refer as the 
completion of computation problem appears: we are never sure to find the computation com- 
pleted, however carefully we choose the instant of time at which to measure the clock. Not 
only: the clocking excitation goes back and forth along the program line, thus doing and 
undoing the computation. On the other side, the entanglement between the register and the 
cursor makes it possible to concentrate all the quantum uncertainty about the state of the 
computation on the position of the clocking particle [l5| . 

Given the initial condition for the program line \ Q = xq) , with xq G {1, . . . , s}, by taking 
the limit of a "long" computation (s oo) we obtain 



lim c(i, x; s) 



lim exp 

k=i 



iXt cos 



kiT 
s + l 



sm 



/ kTTXo \ 

\s + l 



sm 



kiTX 

s+l 



— / exp(iAcosy) sin(yxo) sin(xy)dy 

TT 



(3.12) 
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Jy{t) being Bessel functions of the first type [2a] and A = 1. In particular, for xq = 1 we get 



i"-^(J,_i(t) + J,+i(t)) = i 



(3.13) 



It is worth noticing that the solution in the s — > +oo limit gives a very good approximation 
of the behavior of the system even in the case of a finite chain as long as we consider periods 
of time sufficiently short not to include reflections of the wave packet associated to the cursor 
on the boundaries (see figure 13.20 . 



P(10,t) 
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Figure 3.2: solid red line: |c(i, 10; 20)p; dashed blue (^Jio(^))^- The approximation via Bessel func- 
tion is precise until the wave packet gets reflected on the boudary of the approximated finite chain. 



3.2 Relaxing the upper bound 

The main purpose of this Section is to give examples of the behavior recalled in the assertions 



and (ii) made at the beginning of this chapter. 



This we do with the help of the following Hamiltonian: 

s-l 

i=i 

+p_r+(s + l)r_(s) + 

s+S-l 

+ r+(j + lK(i)+/i.c.). (3.14) 

j=s+l 

Uj, 1 < J < s — 1 being unitary operators acting on the register. With respect to the 
Hamiltonian ()1.2p . we have introduced an additional control q-bit p = {pi,p2,P3) in the term 
r+(s + l)p_r_(s); this is an example of a conditional jump in the quantum walk performed 
by the cursor: it acts non trivially only in the eigenspace belonging to the eigenvalue +1 of 
P3, enabling the transition \Q = s) ^ \Q = s + 1). If this transition is enabled, then the 
cursor can visit the additional telomeric sites s + 1, . . . , s + 5, else it gets reflected back. 



3.2. RELAXING THE UPPER BOUND 



27 



Figures 13.31 and 13.41 give examples of the behaviour of the probabiUty 



s+5 



(3.15) 



j=s 



of finding the register in the state A = Us~i U2 ■ Ui\ R{1) ), under two different initial 

conditions, which determine two different forms of the amplitudes 7. 

Figure 13.31 corresponds to the initial condition | ^1 ) = \p3 = — 1)(8>|(5 = 1)^ the motion of 
the cursor remains confined to the sites 1, . . . , s, as it is 7(t, k) = c{t, k; s) if 1 < k < s, 
otherwise. The probability -P(s<Q)(i) of finding the computation completed satisfies in this 
case the inequality: [2^ 

(3.16) 



Pis<Q)it) < — 



S3 



Figure 13.41 corresponds to the initial condition: | ipi 



\P3 



+1)0 \Q = 1), leading to 
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Figure 3.3: = |p3 = -1) ®\Q = 1); s ^ 20. 



7(t, k) = c{t, fe; s + (5) for 1 < < s + (5. For t just below s + 26 the probability -^^Q)(t) of 
finding the computation completed is close to the much less severe upper bound: [2a| 



1-^f. 



1 



1 + 26/s 



(3.17) 
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Figure 3.4: | V-i ) = |p3 = +1) ® |Q 1); s = 20; S = 10. 
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3.3 The Quantum END Instruction 




t 



Figure 3.5: The solid line represents the probability of finding the cursor in the telomeric chain using 
a IT -pulse applied at time tQ = s + 25. The dashed lines correspond to Fias \3. 3\ and \3.4\ 

The abrupt collapse of P(^g^Q-^{t) at time t « s + 25, evident from Fig. 13.41 corresponds 
to the following fact: traveling with average speed close to 1, at time t ^ s + 26 the cursor 
"returns down the active part of the program line", thus, in effect, undoing the calculation. 
Bringing the computation to an END, and storing the result is not completely trivial in 
the case examined here of a reversible quantum clocking mechanism: "Surely a computer 
has eventually to be in interaction with the external world , both for putting data in and for 
taking it out [3]". 

A simple model of such interaction is suggested by inspection of the Hamiltonian (j3.14p : 
starting from the initial condition I/J3 = +1), the transition \Q = s) ^ \Q = s + 1) is enabled 
by the control term r+(s + l)/9_r_(s) which, simultaneously, determines the transition |p3 = 
+1) ^ I/53 = -1). 

The transition \Q = s+1) — > \Q = s), enabled by the hermitian conjugate term r+(s)/9-(_r_(s+ 
1), will be therefore inhibited if, by external means, we enforce the transition I/93 = —1) — > 
\p3 = +1) at a time, close to to = s + 25, when most of the probability mass is in the region 
s,. . . ,s + S. 

Figures 13.51 (where Figs. 13.31 and 13.41 are also reproduced for comparison purpose) presents 
the effect of the addition to (I3.14p of the time dependent perturbation 

h{t) = B{t) ■ pi (3.18) 

where the "magnetic field" B{t) is non vanishing only in a unit time interval around to, in 
which it takes the value vr. 

With a probability depending only on the ratio 5/s (see ()3.17p ) between the lengths of the 
active part and the telomeric part of the program line, the 7r-pulse (j3.18p definitively prevents 
the cursor from undoing the computation. 

The idea of a n-pulse trap just presented works only if the control q-bit is initialised in the 
\p3 = +1) state. It is immediate to convince oneself that the following double trap Hamiltonian 
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does not suffer from the above limitation: 
H = h{t) + 

s-l s+<5-l s+2(5-l 

-{^U,T+{j + l)T^{j)+ T+{j + l)T^{j)+ r+(i + l)T_(j) + 

j=l i=s+l j=s+(5+l 

+p_ T+(s + 1)t_(s) +P+ T+(s + (5 + l)r_(s) + /I.C.). (3.19) 

With any initial condition for the control q-bit, under the action of the above Hamiltonian, 
the 1/93 = +1) component of the state gets definitively trapped in the first telomeric region 
{s + 1, . . . , s + 5}, the |/>3 = —1) component in the second one {s + (5 + 1, . . . , s + 25}. 

As a final remark of this section, we observe that, acting in effect as a Stern-Gerlach apparatus 
providing space separation between two different spin states, the term 

switch = {p-T+{s + l)r_(s) + p+T+{s + 6 + l)r_(s)) + h.c. (3.20) 

can be used also to model the preparation ("putting the data in") of a register qubit in a 
given spin state. 



3.4 Synchronizing the computational paths 

When defining the C^NOT circuit (section 12. 3p we made massive use of delay lines, that is 
edges of the clocking device "during which" the state of the register does not change. As we 
said in there, this allows for every computational path to be of the same length q in turn, 
this makes it possible the interference between different computational paths corresponding 
to different conditions of the controlling qubits to happen. 

As an example of the role of synchronization in our interacting quantum walks, let us consider 
the following Hamiltonians. 



Hi = 


SWITCH{1,6, 


0-3(1^), 1 




^-(1)|2)(1| + 




+ 


asiiy)\3){2\ + 




+ 


(T+(l)|6)(3| + 




+ 


^+(1)|4)(1| + 




+ 


|5)(4| + 




+ 


a_(l)|6)(5| + 


H.c, 


and 






H2 


= ^-(1)|2)(1 


1 + 




+ ^^3(^^)1 3 )( 2 1 


+ 




+ ^+(1)|5)(3 


1 + 




+ ^+(1)|4)(1 


1 + 




+ ^-(1)|5)(4 


1 + H.c. 



(3.21) 



(3.22) 



^The same principle is applied in electric circuits when two or more signals are required to arrive at the 
same time through paths of different length. 
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(r3(2) 



R(1) 




Figure 3.6: Synchronized circuit 




Figure 3.7: Non synchronized circuit 

corresponding to figure 13.71 

The difference between the two Hamiltonians consists in the lack of the delay line 4 — > 5 in 
(13:22]) . 

We set the initial state of the register to 



i?(l)) = I ai(l) = +1,(73(2) = -1 



(3.23) 



and the initial condition for the cursor to | C(l) ). 
The initial state of the machine is 



4>,) = \R{l))\C{l)). 



(3.24) 



We solve the Cauchy problems 



^'t))=H^\<p{t)) 



dt 



m ) = I '/'I 



(3.25) 
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Figure 3.8: Synchronous case. The probability of finding, upon measurement, the cursor at position 6 
and the register in the state \ o'i(l) = +1,(73(2) = —1). 



U\<t^\t))=H,\c^'{t)) 
ll'A'(0)> = |<Ai) 

The choice of the initial state of the machine is such that both branches of the switch are vis- 
ited with the same probabihty amphtude. In the circuit of figure [3^ along the upper branch 
the operator (73(2) gives a ivr phase factor when acting on the register; the lower branch does 
nothing. Since the upper and lower computational paths have the same length, the clocking 
signal are synchronized and interfere destructively at site 6 giving the result of figure 13.81 
As figure [3^91 shows, the situation is completely different in the case of the circuit of figure [3^71 
In fact, the computational path along the upper branch is one step longer than the compu- 
tational path along the lower branch. Due to the lack of synchronization the interference 
pattern of the clocking signals at site 6 is completely different from the one of the synchro- 
nized circuit. 

The interference between different computational paths plays a key role in quantum algo- 
rithmics which, in fact, is essentially based on interference and on entanglement [10|]. As 
anticipated in section 12.31 this architectural choice makes it possible to define the basis of 
logical successors of a given initial state of the machine in the same way adopted in the case 
of the linear program line. As a matter of fact, the Peres basis for the interacting quantum 
walk of Feynman is a generalization of the projection mechanism used in [ill to study the 
motion of a quantum walker on a graph G„ such as the one of figure 13. Ill Similar results 
hold in the case, studied in [22], in which, because of conditional jumps in the program line 
(such as the ones needed in the iteration of quantum subroutines or the implementation of 



the C^NOT), the cursor performs, in effect, a continuous time quantum walk [2ll . l22l | on a 
planar graph. Indeed, if we traced out the register degrees of freedom we could take the 
basis 

|co/(i)) = -^ E l^ = «) (3.27) 



aG slice j 



where Nj is the number of program line sites in the j-th slice; the vectors | col{j) ) are hence 
uniform superposition of the basis states of TCcursor belonging to the j-th column. With this 
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P..(t)(Q=6 A(r3(1)=-1) 




Figure 3.9: Asynchronous case. The probability of measuring Q — 6 and the register in the state 
I ai(l) = 1,^3(2) = -l>. 




cr3(v) 



R(1) 




Figure 3.10: 
or levels. 



The same SWITCH of figure 



The vertical dashed lines divide the circuit 



slices, 



notation, the motion on the SWITCH circuit of figure [3TT0] of the clocking excitation is exactly 
the same as the one of the traveling excitation of 21 1 . 



Summary 



The qualitative statement "the longer is the computation the smaller is the probability of get- 



ting the cursor in the END position" has been made quantitative in 26|]: Popt{s) = 0{s~ 
By means of a telomeric chain of length (5, that is extra space resources, the latter bound can 
be raised to 0{\/5/ s). We proposed a measurement scheme which has the same probability 
of success; it exploits the possibility of controlling the trajectory of the clocking excitation by 
means of the state of a register qubit: we wait for most of the probability mass of the cursor 
to be in the telomeric region and than prevent the excitation to get out of the telomeric chain 
by disabling the transition s ^ s + 1. 
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Figure 3.11: The graph G4. 



In Feynman's scheme the functioning of a quantum algorithm gets a "topological" char- 
acterization: different computational paths are simultaneously visited, each with a certain 
amplitude depending on the initial state of the register; at a certain point (in space) the 
signals coming from alternative routes interfere with each other. To get the right interference 
pattern it is necessary that the clocking signals are synchronized with each other. 
The synchronization of different computational paths, on the other hand, makes it possible 
to realize the projection mechanism which, in turn, allows for the motion of the cursor on 
planar graphs to be described as if the clocking agent were traveling on a linear chain. 
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Chapter 4 

Speed, entropy and energy 



In this chapter we study the random variable time of flight speed of the clocking agent on a 
hnear chain. A quantitative assessment on the spreading rate of the wave packet associated 
to the clocking excitation makes it possible to define a quasi-optimal preparation scheme 



and post measurement strategy 30[|. We show, moreover, that the intrinsic uncertainty on 
the position of the cursor induces decoherence on the state of the register and that this 
decoherence is of Lindblad type [sij . 

Throughout this chapter we will consider Hamiltonians of the form 

A 

^ = - 2 E ^-1 + 1) )^ C'(x) I + U-^\ C{x) )( C(x + 1) I, (4.1) 

x=l 

where only Ui,U2, ■ ■ ■ ,Un-i are assigned by the algorithm we are interested in whereas 
Un, . . . ,Us-i are to be assigned as a part of the description of the clocking mechanism. 
For instance, we have seen in section [3^ the case in which N — s = 5, and Un, ■ ■ ■ , Us-i = Ir, 
the identity in TLregister , and verified the role of the cursor sites A^, . . . , s as a storage mech- 
anism of the desired output state | R{N) ). 



An alternative point of view was taken in some of the numerical examples of |25l |. motivated 
by Grover's algorithm: one may suppose all of the Ux to coincide, in such a context, with 
Grover's estimation ■ oracle step G, and study the effect of applying G more than the optimal 
number A^ — 1 of times. 

Most of our numerical examples will refer in fact to the following particular instance (Toy 
model) : 

A 

Ht = --Y^ e-^t'^a \c{x + l)){ C{x) I + e^t'^a | ^{x) ){C{x + l)\. (4.2) 

^ x=l 

which, with the right choice of the parameter a, models the search of the "needle in the 
haystack" of Grover algorithm 



4.1 Speed of computation 

In Chapter [T] we have seen that a Feynman machine initialized (at time t = 0) in the state 

I R{1) ) ® I C{1) ) = I Ml ) (4.3) 
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evolves, under the Hamiltonian ()4.ip . into 

s 

I Mi(t) ) = e-*^*| Ml ) = ^c(t,2;;s)| R{x) ) | C7(a 



(4.4) 



x=l 



where the register states | R{x) ) are defined as in ()1.15p . 

The observable Q/t acquires thus the meaning of number of primitives per unit time applied 
to the initial condition | R{1) ) in the time interval (0,t). In order to study the behavior 
over long intervals of time (t — > +oo) of this observable in the case of a long computation 
(s — > +oo) it is expedient to study its characteristic function 



(l)s,tiz) = {Mi{t) 



exp \ iz^\ I Mi{t) ) = |c(t, x; s)|^ exp {iz— 

^ ^ x=l 



(4.5) 



namely the Fourier transform of its probability distribution. 

Theorem 4.1.1. ^ 

lim lim (l)s,t{z) = / 



(4.6) 



Proof. The large s behavior is easily studied by inserting the explicit integral representation 
of the s +00 hmit of (EJl) into KB 



lim c{t, X] s) 



Jt cos(p) 



sin(p) sm{xp)dp. 



We thus obtain 



lim > 



e t 



x=l 



- I e 

vr 



it cos (p) 



sin(p) sin(xp) dp 



(4.7) 



(4. 



The t — > +00 limit is studied by substituting the sum over v = x/t, step 1/t, appearing in 
(14. 5 p with an integral 



lim 



2Vt r 



vr 



Jtcos{p) 



s'm{p) s'm(vtp) dp 



dv. 



(4.9) 



We evaluate the leading contributions to (j4.9p by a standard stationary phase argument. 
First of all we observe that 



Jtcos{p) 



sin(p) sm{vtp) 



ivtp I -ivtp 
^ifcos(p) I ' 

2i 



pit(cos{p)+tip) _j_ git(cos(p)— vp) 



sin(p) 



2i 



■ sin(p). 



(4.10) 



Since the equation 



dp 



(cos(p) + vp) = 



^The same result (see also [2^) has been obtained by N. Konno [33] in a different context: there the author 
investigates the propagation of an excitation through a homogeneous tree, whereas here we deal with a linear 
chain. The obvious explanation of the coincidence of the results in two such different contexts comes from the 
projection mechanism explained in section \3A\ 
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admits two distinguished solutions 



pi = arcsin{y) (4-11) 

P2 = TT — arcsin{v) (4-12) 

whereas the equation 

— (cos(p) — vp) = 
dp 

has no solution for < v < 1, the contributions to the limit 

(riT ^it{cos{p)+vp) 2^ /"I" git(cos(p)-vp) \ 
/ {p)dp + I sin(p) dp (4-13) 
Jo sm Jq 2i J 

come only from the first term; since these contributions come from the region around pi and 
P2, the following equality holds 



/■TT At{cOs(p) + Vp) /■ + 00 

lim / dp= lim y (^ii^°<Pj)+m) / e-^^^^'^P^^^-P^) dp. (4.14) 



rTT pit{cos{p)+vp) 

■dp = lim 

ie{i,2} 

We recall that 



1 



mil ( — y= e 5^ ) = 5p^{p) (4.15) 



where Sp. is the Dirac-(5 function concentrated around pj . 
The right hand side of equation (|4.14p is then equivalent to 



lim y e^*(™"fe)+^P^) J . (4.16) 

*""+~jg7t2} ^ itcosipj) 

By inserting the last expression in (14. 9p we get 

^izv 4(t;^ ^ ^.^ /'^ ^i^,; C0S((C0S(P1) - C0S(P2)) + Vjpi - P2))) ^^.^ 

vrVl — ^^+°°Jq sJ COS (pi) cos(p2) 

Since the equation 

— (cos(pi) - cos(p2)) + v{px - P2) = (4.18) 
dv 

has no solution for < f < 1, the second term vanishes and we get the result. □ 

As convergence in the sense of characteric functions implies convergence in the sense of 
cumulative distribution functions (^convergence in law), we conclude that a "long" computa- 
tion starting from the initial condition (14.31) proceeds "in the long run" at a rate of V{Mi) 
steps per unit time (the unit of time having been set so that A = 2 ), V{Mi) being the 
random variable defined by having as its characteristic function the right hand side of ()4.6p : 
equivalently stated it has probability density 

4t;2 

fv{M,){v) = 7(0,1) (^) f (^-19) 



vrVT" 



V 
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Here and in what follows we denote by /(a,fe)tlie indicator function of an interval (a, b): 

, , . f 1 if X G (a, 6) r,r>\ 

'iaM-) = \ otherwise. ^''''^ 



The mean value 



and the variance 



E{V{M,)) = j\ fv(M,){v)dv = ^ (4.21) 



var{V{M,)) = E {{V{M,)f) - E {V{Ivh)f = ^ " (|^) (4-22) 
are then easy to compute from (I4.19|) . 

It is worth extending our analysis to more general initial conditions; for any positive integer 
xo, a state such as 

I = I R{xo) ) ® I C(xo) ) (4.23) 

having at a certain instant the cursor in xq, evolves under the action of the Hamiltonian (j4.ip 



as 



I M^At) ) = e"'""*! M(Xo) ) = ^c,o(t,x,s)| R{x) ) ® \ C{x) ). (4.24) 

x=l 

With the same techniques used in the proof of Theorem 14.1.11 we can prove that 

Theorem 4.1.2. The random variable time of flight speed V{Mxq) has a characteristic func- 
tion 

/xi7V<3i,^ ,^^ 4, s\\i{xq arcsiuiv)'^)'^ , . 

lim lim ( M^o(t) \e''T\ M^,{t) ) = ^ . ^ -■ (4.25) 

The cumulative distribution function of V{Mxg) is consequently 

FviM.,){v) ^ Prob{V{M,J<v)= (4.26) 
, , /2arcsin(f) sin(2xo arcsin(w)) \ ^ , , 

= ^(0.1) — — — + ^(l,+oo) 2; 

corresponding to an expectation value 

E{V{M,,)) = ^ (4.27) 
4tt — tt/xq 

Comparison between (|4.2ip and (|4.27|) shows the effect of a measurement of Q. If, at a given t, 
Q is measured and the result xq is found, then the state ()4.4p . into which the initial condition 
(|4.3p has evolved, collapses into the state (|4.23p . From this moment on the computation 
proceeds at the mean rate (I4.27p : for large values of t, reading the clock is likely to reduce 
the speed of further computation by a factor 3/4 (without, because of (j3.6p . altering its 
correctness). 
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4.2 Entropy 

Motivated by the experience gained under the particular initial conditions (j4.3p and (j4.23p 
we define, for any (unentangled) initial condition of the form (for fixed e > 1) 

€ 

\R;^l^o) = \R)^Y.'f'o{x)\C{x)), (4.28) 

x=l 



the "time-of- flight speed" 3j] of computation in the state V'o as the random variable V{ip[ 



0) 



having characteristic function 

cf>v(^bo)i^) = lim lim ( R^tPo |e'*^exp (iz^) e"**^! R■,^l^o )• (4.29) 

^ ' t^+oo s^+oo \ t J 

The above limit is easily shown to exist by the techniques outlined in the previous section; it 
corresponds to the probability density 

|^(arcsin(7;))p + [^(tt - arcsin(t;))p 

fv(^Po) [V) = 1(0,1) (v) (4.30) 



where 



"^{p) = ^J^^sm{px)Mx) (4-31) 

x=l 

is the sine transform of the initial state of the cursor. 

The observable Q retains in this context the meaning of relational time [sS] in the sense that, 
given that at any parameter time t the cursor is found at x, it is then certain that the register 
is found in the state Ux-i • . . . • U2 • Ui \ R) . 

In reading the output at any time t, namely in the measurement of any, however carefully 
chosen, observable of the register, there is an intrinsic uncertainty corresponding to the un- 
certainty about how far the computation has proceeded. The fact that Q/t has a non trivial 
limit in law means that the leading term of the variance of Q is proportional to t'^ and there- 
fore that the uncertainty increases with t. This section is devoted to the examination of an 
example in which the notion of "the most careful choice" of the observable to read on the 
register can be made precise and shown to be pertinent to the algorithm considered. 
We consider for the moment the initial condition | Mi ) given in (14. 3p and its time evolution 
I Mi{t) ) described in (j4.4p . More general initial conditions of the form (|4.28p will be examined 
in the next section. 
Call 

Pm{t) = \Mi{t)){Mi{t)\ (4.32) 
the density matrix of the machine at time t. 

By taking the partial trace Tr-j-i^^^^^^ {pm{t)) with respect to the cursor degrees of freedom, 
we get the density matrix pr{t) of the register: 

s 

p^{t) = ^|c(t,x;s)p| R{x) ){R{x) \. (4.33) 

x=l 

Call Xj{t) the positive eigenvalues of Pr{t) and \ bj{t)) the corresponding eigenstates. A simple 
computation, amounting to the Schmidt decomposition [s^ of the state (14. 4p . shows, then. 
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that the density matrix of the cursor is given by 

Pc{t) = Y.XMdj(.t)){d,{t)\ (4.34) 
j 

where 

I d,it) ) = -A^j2<t,x;s){ b,{t) I R{x) ) I C{x) ). (4.35) 

x=i 

Because of (I4.34p and of the orthonormahty of the states | dj{t) ), the von Neumann entropy 
of the register and also of the cursor is then given by 

S{Pc{t)) = A,(t)lnA,(t) = S{pr{t)). (4.36) 

j 

We observe that, as (I4.33P shows, the von Neumann entropy of each subsystem does depend 
on the algorithm being performed. It is, indeed, only under the hypothesis, nowhere made 
above, that the states | R{x) ) are orthonormal that ()4.33p is the spectral decomposition of 
Pr{t) (the von Neumann entropy becoming in this case equal to the Shannon entropy of the 
distribution of Q). 

We focus our attention, in what follows, on our Toy model ()4.2p . in which the register is a 
single spin 1/2 system. We indicate by 61,62,63 the versors of the three coordinate axes to 
which the components a = (o"i, o"2, (T3) of such a spin are referred. 
In the basis | CJ3 = ±1 ), the density operator Pr{t) will be represented by the matrix 

p(t)--( ^ + '3^') S,{t) - Z S2{t) \ 

where 

Sj{t) = Tr {pr{t) ■ Gj) , J = 1, 2, 3. (4.38) 

Equivalently stated, the Block representative of the state Pr{t) is given by the three-dimensional 
real vector 

s 

s{t) = ^\c{t,x;s)\^ {R{x) \a\ R{x) ). (4.39) 

x=l 

We shall assume, in what follows, that the initial state of the cursor is | C(l) ) and that the 
initial state of the register is of the form 



I R{1) ) = cos l^-J I ^73 = +1 ) + Sin (^-J I (73 = -1 ) (4.40) 

namely the eigenstate belonging to the eigenvalue +1 of n(l) - a, with 

n(l) = 61 sin 6* + 63 cos 6*. (4.41) 

We wish to remark that the above example captures the geometric aspects not only of such 
simple computational tasks as NOT or V NOT (viewed as rotations of an angle vr or 7r/2 
respectively, decomposed into smaller steps of amplitude a) but also of Grover's quantum 
search js^]- If, indeed, the positive integer // is the length of the marked binary word to be 
retrieved, and we set 

x{p) = arcsin(2"2 ) (4.42) 
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and 

d = TT-2 xip) (4.43) 

then the state ()4.40p correctly describes the initial state | ) of the quantum search as having a 
component 2~^/^ in the direction of the target state, here indicated by | a; ) = | (J3 = +1 ) , and 
a component ^/Y^^2~^ in the direction of the flat superposition, here indicated by | (73 = — 1 ), 
of the 2^^ — 1 basis vectors orthogonal to the target state. In this notations, if 

a = -4x(/u), (4.44) 

then the unitary transformation exp(— zacJ2/2) corresponds to the product B- Aoi the oracle 
step 

A = Ir-2\uj){uj\ (4.45) 

and the estimation step 

B = 2\l){l\- Ir. (4.46) 



Following the beautifully pedagogical approach of Jozsa [37|, we observe that the operator 
A corresponds to a reflection in the hyperplane orthogonal to | cj ) and B corresponds to a 
reflection in the hyperplane orthogonal to | i ), with a minus sign. 

Lemma 4.2.1. If \x) ^^1/ state in firegister , then = — 2| ■0 )( -0 | preserves the 
subspace S of TCregister Spanned by \ x) and 

Proof. Geometrically, S and the mirror hyperplane are orthogonal to each other (in the sense 
that the orthogonal complement of either subspace is contained in the other subspace) so the 
reflection preserves S. Alternatively in terms of algebra, ()4.46p shows that B takes \ t) to 
— and, for any | a; ), it adds a multiple of | w ) to| t ). Hence any linear combination is 
mapped to a linear combination of the same two states. □ 



Since Grover's algorithm consists of the iterated application of two reflections (a reflection 
in the hyperplane orthogonal to the target state | uj ) and a reflection in the hyperplane 
orthogonal to the initial state | i )), it preserves the space S spanned hy \ u) and \ l). 
Now we may introduce a basis {| ei ), | 62 )} into S such that | u ) and | t ), up to an overall 
phase, have real coordinates. Indeed choose | ei ) = | z. ) so | i ) has coordinates (1,0). Then 
e^^\ u>) = o| ei ) + 6| 62 ) where | 62 ), orthogonal to | ei ), has still an overall phase freedom. 
We can thus choose ^ to make a real and the phase of | 62 ) to make b real. In this basis, then, 
since | lj ) and | l ) have real coordinates, the operators A and B when acting on S, are also 
described by real 2 by 2 matrices: they are just the real 2 dimensional reflections in the lines 
perpendicular | uj ) and \ l). Finally we have the following 

Theorem 4.2.2. Let Mi and M2 be two mirror lines in the Euclidean plane M intersecting at 
a point O and let a be the angle in the plane from Mi to M2 . Then the operation of reflection 
in Ml followed by reflection in M2 is just rotation by angle 2a about the point O. 

It is having in mind the connection with Grover's algorithm that, for the sake of deflnite- 
ness, in the examples that follow we are going to consider the one-parameter family of models, 
parametrized by the positive integers fx, corresponding to the choice ()4.43p and (j4.44p of the 
parameters 9 and a and to the choice s = 2^ + 1 of the number of cursor sites, corresponding 
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S3 




Figure ^.1: A parametric plot of (si(t), S3(t)) for Q <t < s, A = 1. The choice ^ = 7, x = 
arcsin(l/2''/^), s = 2^ + 1, a — — 4x, = tt — 2x of the parameters is motivated by the connec- 
tion with Grover's algorithm. Only the initial state lies on the unit circumference, the locus of pure 
states. 



to the possibility of performing up to an exhaustive search. 
In the example defined by the above conditions it is 

( R{x) I a\ R{x) ) = sin {9 + {x - l)a) + cos {9 + {x - l)a) 63 (4.47) 

and, therefore, 

s 

s{t) = ^ \c{t, x; s)|^ (sin {9 + {x - l)a) + cos {9 + (x - l)a) 63) . (4.48) 

x=l 

Figure 1 presents, inscribed in the unit circle, a parametric plot of (si(t), S3(t)) under the 
above assumptions . It is convenient to describe the Bloch vector s{t) = si{t) Ci + s^^t) 63 in 
polar coordinates as 



si{t) = r(t)sin7(t), S3{t) = r(t)cos7(t). 



(4.49) 
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A very simple approximate representation of s{t) becomes then possible: 

s 

r(t)e*^W = ^|c(t,x;s)pexp(i(0 + (x-l)a)) = 

x=l 

s 

= exp(i(0 — a)) \c{t, x; s)p exp(iax) 

x=l 

^ exp{i{e -a))E{eyip{ia\tV{Mi{0)))) (4.50) 

The last step, legitimate for 1 << A t < s, requires only the explicit computation of the 
characteristic function corresponding to the probability density (I4.19p . which leads to 

^^tyy{t) ^ ^^^P^^(^^~"^) ((Ji(r) -TJ2{T)+i{THo{T) - //i(T)))) (4.51) 

where Jk and Hf^ are, respectively, Bessel functions and Struve functions and T = aXt. 
The time evolution of the register subsystem is summarized by the Lindblad equation 3^, 31] 



dprit) i d'yit) r , ldlnr(t) , , , ,^ ^. 

The commutator term [a2, Prit)] describes the Hamiltonian part of the dynamics (after all we 
are considering a rotation about the X2 axis); the double commutator [(J2, [cj2, /Or(i)]] describes, 
in much the same sense as equation 2.8 of the decohering effect of this rotation being 
administered by the cursor in discrete steps at random times. 
The eigenvalues of Prit) can be written as 

Ai(t) = i(l + r{t)), X2{t) = i(l - r{t)). (4.53) 

The von Neumann entropy S (pr(i)) is therefore 

S = _l±ZW,„i±^_iz^„iz^, (4.54) 

An example of its behaviour is shown in figure 14.21 The eigenvectors corresponding to the 
eigenvalues (14.530 are, respectively 

I hit) ) = ( ^'f^'jiii ) , I b2{t) )=(- ^-(Ti^y^) \ _ 



sin(7(t)/2) ; ' ' ^ ' V cos(7(t)/2) 

It is to be stressed that, at each time t, the projector | bi{t) ){ bi{t) \ is, among the projectors 
on the state space of the register, the one having in the state Pr{t) the greatest probability of 
assuming, under measurement, the value 1. Thus, the most careful choice (the one affected 
by minimum uncertainty) of the observable to read on the register at time t is the projector 
I bi{t) ){bi{t) \. In the case of Grover's algorithm one must measure the projector | u; )( a; | = 
I CJ3 = +1 )( CJ3 = +1 I (and one easily can, because of the kickback mechanism analyzed, 
for instance, in [25>] ) and has the freedom of choosing the time r at which to perform the 
measurement. The best choice is therefore such that | &i(t) ) = | cj ) (in our notational setting, 
r is the time at which the helix of figure WA\ crosses for the first time the positive S3 axis). 
In spite of the fact of being now in the most favorable setting, one has, nevertheless, a deficit 
1 — Ai(r) in the probability of finding the target state. 
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S(Pr(t)) 




Figure 4-3: The same model as in figure and figure < t < 1.2 s. The thin solid line is a 
graph of Tr{pr(t) ■ {Ir + (T3)/2), the probability of observing the target state \ uj ) = \ ^ +1 ) in the 
example of Grover's algorithm. The dashed line is a graph of Tr{pr{t) ■ [Ir — i73)/2), the probability 
of observing the "undesired" output \ = — 1 )• The upper and lower bounds on the probability of 
observing the target state are represented by the thick solid lines Xi{t) and X2{t). 



As figure [^31 shows, there are successive instants of time at which the probability of successful 
retrieval has a local maximum (a remnant of the periodic nature of Grover's algorithm when 
applied by an outside macroscopic agent) but the heights of these successive maxima form a 
sequence having a decreasing trend. 

Further insight into our toy model is gained by examining the t dependence of E{Q{t)) = 
{ Mi{t) \Q\ Mi{t) ) and of the angle of polarization j{t). The example of figure IT^ a suggests 
that the mean value of speed derived from asymptotic considerations correctly describes the 
average behavior of the "clocking" subsystem also for finite values of < t < s. As figure HT^ b 
shows, the "clocked" subsystem system cr is, in turn, driven, on the average, into uniform 
rotational motion. 

Our model is so simple that we can explicitly study how the above semiclassical picture 
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(in which the time parameter t acquires operational meaning from its hnear relation with 
mean values of configurational observables of clocking and clocked subsystem) is distorted 
by a measurement performed on either subsystem. The observations made at the end of the 
previous section about the effect of reading the clock can indeed be complemented by the 
examination of the effect of reading the register. 

Suppose that the observable cr^ has been measured at time r and the result +1 has been found: 
the Bloch diagram of figure I4.5l a shows then that the evolution of the register proceeds in 
much the same way as in the undisturbed situation of figure 14.11 (with the only obvious 
difference that the post-measurement initial condition | 6i(t) ) lies on the unit circumference). 
If, instead, the result —1 has been found (figure H31 bl the post-measurement evolution of the 
register is completely different from the unperturbed one. 

We conclude this section with an example of the insight that the time evolution of S{pr{t)) 
can give on the algorithm Ui,U2, . . . , Us-i being performed by the machine. Suppose of using, 
instead of the assignment (14. 2p of the primitive steps, Ui = U2 = . . . = Ug-i = exp(— iQ:a2/2), 
the alternative assignment 

{A for odd x , , 

Ux = \ „ f (4.56) 
[ B for even x 

where A and B are given by (j4.45l) and (j4.46p . Figure gives, for this example, a full 
account of the diffusive character [23|] of Grover's quantum search: the first maximum of 
the probability of finding the target state (figure 14.61 a) is reached in correspondence of the 
first local minimum of entropy (figure 1321b): that the search has gone, before this instant, 
through a local maximum of entropy is shown with particular evidence by the Bloch diagram 
of figure [Ol c. 



4.3 Energy 

With reference, for definiteness, to the example of figure 3, call r the instant of time at which 
the probability Tr{pr(t) | a; )( a; |) reaches its first and absolute maximum. We recall that, 
in the above example, the target state | a; ) is taken to be the "up" state | 173 = -|-1 ) of the 
register. 

The whole point of the analysis of the previous section is that Ai(r) is strictly smaller than 1. 
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0.2 0.4 0.6 0.8 1 1.2 

(c) 

Figure ^.5: Same parameters as in figure [^T7| ' measurement of the observable at time r. Frames 
(a) and (b) represent the evolution for t < t < At of the Bloch vector when measurement returns +1 
and —1 respectively. Frame (c) represents the cumulative distribution functions of the speed V when 
the results +1 (solid thick line) and —1 (dashed line) have respectively been found; the solid thin line 
represents the c.d.f. of V in case of no measurement. 



This amounts, in turn, to a deficit 1 — Ai(t) in the probabihty of finding the target state. 
This deficit is not, in itself, a strong limitation in a quantum search algorithm, because we 
can in principle identify the right target through a majority vote among a "gas" of a large 
number N of machines. The trouble is that if we want to use the same machines once more, 
we need to purify the "gas" of registers from the fraction AoIt) of them which have collapsed 
into the wrong state: standard thermodynamic reasoning [411] shows then that this requires 
the removal from the gas, supposing a heat reservoir at temperature T is available, of an 
amount of heat of NkBTS{pr{T)), ks being Boltzmann's constant. 

We wish, in this section, to supplement the above considerations with an explicit description 
of the post-measurement state of the machine, showing, in particular, the effect onto the clock 
of the act of reading the register 4^] . 

Suppose that at the optimally chosen instant r, at which it is 7(t) = 0, while the machine is 
in the state | M(r) ), a measurement of the projector [I^ + C3)/2 is performed. 
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1 r 




100 s 150 200 

(a) 



100 s 150 200 
(b) 




(c) 

Figure 4^.6: Same parameters as in the previous figures; Ux as in \4.5b\ l. The Bloch diagram refers to 
the time interval (0, t) needed to reach the first maximum in probability. We observe that the state of 
the register passes through the maximally mixed state before cooling down toward the target state. 



1 

as = +1 ) V c(r, X- s) cos{{0 + (x - l)a)/2)\ C{x) 



If the measurement gives the result 1, then the state | M(r) ) collapses to 
|Mi(r)) 

If, instead, the measurement gives the result 0, then the state | M(r) ) collapses to 

\M2{t)) 

Figure 14.7( a) and figure 14.7( b) show the probability distributions 



(4.57) 



0-3 



1 " 

-j== ^ c(r, x; s) sm{{d + {x - l)a)/2)\ C{x) 



Pi{x,t) = |(c(r,x;s)cos((^ + (x-l)a)/2))|VAi(r) 
P2{x,t) = \{c{t, x; s) sm{{e + {x - l)a) /2)) f /X2{t) 



(4.58) 



(4.59) 
(4.60) 



of the observable Q (position of the cursor) in the states | Mi(r) ) and | M2(r) ), respectively. 
Figure Hil]^ a) and figure HiD^b) show the probability distributions of the observable H (energy 
of the machine) in the states |Mi(r)) and |M2(r)), respectively. The two energy distributions 
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Pl(x,T) 

0.3 r 



0.2 



0.1 







15 20 



(b) 

Figure 4-7: Figures (a) and (b) represent, for the same choice of parameters as in figure \4.1\ respectively 
the probabilities Pi{x,t) and P2{x,t) as functions ofx. 
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Figure 4-8: Figures (a) and (b) represent the probability distribution pi[Ek) and p2{Ek) of the energy 
H in the state \ Mi(t) ) and \ A/2(t) ) respectively. 



of figures 5 are easily derived from the fact that the Hamiltonian H defined in (4) has the 
eigenvalues 

Ek = -A cos(i?(A;; s)), k = l,...,s; (4.61) 

each doubly degenerate, an orthonormal basis in the eigenspace belonging to the eigenvalue 
Ek being given, for instance, by the two eigenvectors 

s 

I ^fc; (72 = ±1 ) = I as = ±1 ) ® ^ Vk{x) exp{^ia{x - 1)/2)| C{x) ), (4.62) 



x=l 



where 



Vk{x) = \l sin(x i^ik; s)). (4.63) 
" s + 1 



This leads to the explicit expressions 

PjiEk)= \{Mj{T)\Ek;a2 = v)\\ i = l,2. (4.64) 

7?=±1 

Figures HiH] and [131 show that a collection of identically prepared and independently evolving 
machines becomes in fact, under the operation of reading the register at time r, a mix- 
ture of two distinct "molecular" species, "1" (present in a concentration Ai(t) ), and "2" 
(present in a concentration A2(t)). In each of these two molecular species, the same "atomic" 
constituents have arranged themselves in a different geometrical shape (figures l4.7p . with a 
different orientation of the register spin (equations ()4.57p and ()4.58p ). because of a different 
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energy distribution (figures! 
Comparison with the distribution of H in the pre-measurement state | M(t) ), given in fig- 
ure 14.9^ shows that the presence of the impurities of type "2" is due to unusually intense 
exchanges of energy between the machine and the reading (measurement) apparatus. 



0.01 
0.005 



-0.5 



0.5 



Figure 4-9: In the state \ M{t) ) the probability distribution of H is given by p(Ek) — {vk(x))'^ 



4.4 The role of initial conditions 

An initial condition of the form 

I R{1)- M = \ R{1) ) ® ^o(^)l C{x) ) (4.65) 

x=l 

with ■00 having support in a bounded region = {1, 2, . . . , e} C {1, 2, . . . , s] evolves, under 
dUD as 

e 

e-**^| i?(l); V'o > = Mt, x)\ R{x) ) ® | C{x) ), (4.66) 

x=l 

where ipQ{t,x) solves, with the obvious boundary and initial conditions, the (discretized) free 
Schrodinger equation. The ensuing spreading of the wave packet leads to an increasing trend 
(with the exception of the effects of reflection at time i ~ s evidenced in figures 14. 61 b and 
14. 2p of the von Neumann entropy S{pr{t)) of the state 

s 

Pr{t) = Y,\Mt,x)\^\ R{x) ){R{x) I (4.67) 

x=l 

of the register. This is an undesirable feature because S{pr{t)) gives a lower bound on the 
Shannon entropy of the distribution of any observable of the register, for short on the uncer- 
tainty in any reading of the output. 

The models of the previous section where intended to show the above effect; in this section we 
devote some effort to the goal of decreasing it, by suitable choices of initial condition aimed 
at reducing the spreading of Q in the state 'ipo{t,x). It is sufficient, for this purpose, to study 
only the cursor, evolving under the Hamiltonian 

A 

Ho = --Y,^+{x + l)r_(x) + r+(x)T_(x + 1). (4.68) 

x=l 
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The point is to devise an initial condition tpQ which uses whatever additional finite amount 
Ae = {1, 2, . . . , e} of space resources is available as a launch pad for the cursor in an "efficient" 
way: this means both a high value of the expectation of V{'ipo) and a small value of the 
variance of V^(V'o) (we want the spreading of Q to increase at a low rate for a short time of 
computation). That both goals can be achieved is shown by examining the family of initial 
conditions, given by the eigenstates of a Hamiltonian of the form (14.680 restricted to qubits 
in A^: 

I ) = E VTTT '''' (tTi) ' ^^''^ ), k = l,2,...,e. (4.69) 

x=l ^ ^ 

The probability density of the speed Vk = V{ck) corresponding to each of the above states is 
easily computed from (j4.3Up : 

4 (3 - 2v^ + cos (fl^)) (sin (^)) ' (sin((e + 1) avcsm{v))f 
nVl^ie + 1) {2v^ + cos (f^) - l)' . 

The behavior of Vk is examplified by figure I4.1UI We are taking there, as we will always do 



Fvk(v) 




V 



.28 .56 .78 .92.97 



Figure 4^.10: n = 5, e — 2n~ 1 — The cumulative distribution functions Fv^{v) — ProbiVk < v) 
corresponding to the densities J^. 70^ , for k going from 1 to n. The ticks on the v axes are E{Vi) < 
EiV2) < ... < E(V5). 

in this section for the sake of notational convenience, e to be odd 

e = 2n - 1. (4.71) 

The examples of figure HTQ] clearly show the dispersive nature of the medium (14.680 : they also 
show that increase of the mean value is accompanied by decrease of the variance (as shown 
by the increase in the steepness of the graph as k goes from 1 to n). An obvious choice for 
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the initial state of the cursor emerges from the above example: 

2n-l 



i:/Isin(|.)|C(x)) 

X = l 



C{1) ) - I C{3) ) + I C{5) ) + . ■■ -(-inC(2n-l)) ^ ^^^^^^ 



n 



conforming to the idea of packing the maximum number of wavelengths in the launch pad 
= {1, 2, . . . , e}, and having a, presumably easy to prepare, stationary state of the free XY 
chain localized in A^. 

The random variable Vn = V{cn) has probability density 

_ (sin(2narcsin(z;)))^ 
JVn - ^{o,i)[v) 57:77^ (4.73) 

and, therefore, expectation value 

h=l ^ ' 
h=n+l ^ 

w 1 - — (4.74) 
27rn ^ ' 

The second moment of Vn is explicitly given by 

^(K?) = l-^- (4.75) 

The above considerations lead to the following asymptotic behavior, for large n, of the variance 
of K: 

mr(K) = (4.76) 
47rn 

Equation ()4.74p is a quantitative assessment of the cost in terms of space resources of achieving 
the first requisite of efficiency^ namely high mean speed; similarly, (14.760 gives the cost of 
decreasing the variance of V^. 

Incidentally, as the observable Q has, in the state | c„ ) , expectation value 

E{Qn) = {Cn\Q\Cn)=n (4.77) 

and variance 

r? — \ 

var{Qn) = , (4.78) 

equation (|4.76|) can be read as saying that, in the initial state | c„ ), the position- velocity 
uncertainty product is given by 

var{Qn)var{Vn) ^ ifc^ (4.79) 
izvr 
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E(Q(t)) 




20 40 60 80 100 120 140 

Figure Solid lines: the expectation value of Q in a state |Af„(t)) evolving from an initial condition 

having the cursor in | c„ ); the slope of the initial linear part of the graph is correctly predicted by 
^.74^ . For comparison purposes the dashed lines show ( Mi{t) \Q\ Mi{t) ) as a function of time and 
the corresponding linear fit with slope given by ^.21\ l. 

Figures 14.111 and 14.121 show the relevance of the above asymptotic considerations for the case 
of finite e and finite s for t < s. The effect of the initial condition is most evident if we 
compare the evolution of the state of the register from the initial state | Mi ) = |i?(l))®|C(l)) 
with the evolution starting from 

\Mn) = \R{l))^\Cn). (4.80) 

This is done in figures 14.131 and 14.141 in the same probability-entropy-Bloch format as in 
figure HiHl We examine there two different ways of using an additional amount s — of space, 
of size comparable with the minimum amount required by the algorithm. Figure 14.131 
summarizes the experience developed in [26] on the effect of using all this additional space 
as a telomeric chain or "landing strip" : as long as the cursor stays in this region the register 
remains acted upon by the optimal number of primitives. Figure B. 141 shows the improvement 
obtained by investing part of the additional space as a "launch pad" on which to prepare a 
state in which the spreading of the cursor increases (see figure at a lower rate than when 
starting from position 1. 

Comparison of figures I4.13[ c and 14. 141 c. in particular the improvement of the behavior after 
reflections at site s, shows that the idealized scenario of reversible computation (the cursor, 
"going back and forth" , "does and undoes" the reversible computation) is within reach, with, 
as (I4.76P shows, a polynomial cost in space. We note, in figure 14.151 that we can do much 
better than in figure 14.141 with the same expenditure of space resources, in approximating 
the reversible scenario if, instead of the initial state ()4.72p . we set the cursor in the initial 
state 

x=l 

The state | 7„ ) emerges quite naturally as a three-mode approximation (a linear combination 
of I Cn ) and | c„±i )) of the initial condition that maximizes the mean speed of computation 
for fixed length e of the launch pad. 
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var(Q(t)) 



var(Q(t)) 
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(a) 



(b) 



Figure 4-12: s = 50, n = 5. (a) The variance of Q in a state \ Mn{t) ) as a function oft, compared 
with its best fit of the form const. + t^(4 — tt)/ {Awn), in the time interval (e, s — e) in which boundary 
effects can be neglected, (b ) The variance of Q in the state \ Mi (t) ) , compared with its approximation 
^2(3/4- (8/(37r))2), suggested by UJ^ . 

Summary 

The quantitative characterization of the random variable time of flight speed shows that the, 
classicahy obvious, choice of a sharp initial condition | Mi ) , leads to a fast spreading of the 
wave packet; this in turn, makes the von Neumann entropy of the register grow fast. We 
showed that the entropy of the register is related to the probability of finding the desired 
output state written on the it. We showed, moreover, that we can considerably reduce the 
spreading rate by means of extra space resources, or launch pad, into which we prepare a 
state I 7n which is "the most" localized in momentum representation in a finite region. 
While studying the analytic expression of the entropy of the register in our toy model, pr [t) , 
we gave a simple example of Lindblad dynamics for the state of the system: the clock, which 
administers the computational primitives to the register at random times, acts as a decohering 
environment for the latter. 

In section 14.31 we made quantitative assessments on the energetic perturbation induced by a 
measurement of the system. We observed that, since the Hamiltonian is time independent, 
the only cost of a computation on a Feynman machine is related to the energy required to 
reset the system. What we showed is that it depends on the outcome of our measurement. 
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Figure 4-13: s = 50, <t < 3s, /i — 10, N — [■|2''/^J (the Grover- optimal number of active steps); 
Ui = U2 = ■ ■ ■ ^ Un = exp(— iQ!cr2/2), with a and 9 chosen as in ^-43^ o,nd ^.44^ ; Ux — Ir for 
X > N; initial state \ Ah ) = \ R{1) ) (g) | C(l) ). 
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Figure 4-14-' * — 50, < t < 3s, /i = 10, N = 25, a and 9 as in figure \J.1S\ n — 5, e ^ 2n — 1; 

= Uc+i = . . . = Ue+N-i = cxp(— zacr2/2); Ux = Ir, for l<x<eorx>e + N ; initial state 
|M„) = |i?(l))®|c„). 
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Chapter 5 

A multi-hand quantum clock 



In this chapter we generahze the model of previous chapters by having more than one particle 
traveling on the program line. We will show that, as soon as the operators acting on the 
register do not commute, new effects appear; we provide an example where the excitations 
get confined in bounded regions, an effect resembling Anderson localization 43l. |44 1 . 
The resolution of the eigenvalue problem for what we call the multi-hand quantum clock in 
the non-commutative context of the cursor-register coupling in the general subspace = k 
turns out to be too involved (at least for the author). The numerical examples provided in 
this section are frankly heuristic. Analytic justification must wait for future developments. 



5.1 Number of particles 

In the previous chapter we have provided examples of the benefit of spreading the initial 
wave function of the cursor {N3 = 1) on an initial launch pad instead of, as it would be 
classically "obvious", having it strictly localized at site 1. Equality (j4.76p is, in this context, 
a quantitative assessment of the cost, in term of space resources, of implementing Feynman's 
ballistic mode of computation. 

In this section we abandon, in the same spirit, the classical prejudice of having a single 
clocking excitation, and present a preliminary analysis of the idea of starting the cursor in 
an initial state with > 1. The idea is to follow the motion of a swarm of several clocking 
agents (cursor spins in the "up" state) acting on the register. Stated otherwise, with reference 
for simplicity to the case = 2, we allow the clock to perform a quantum walk on the graph 
having the vertices {xi,X2), with 1 < xi < X2 < s, with edges between nearest neighbors [45}]. 
We recall, mainly in order to establish our notation, a few basic facts [4^ about the XY 
Hamiltonian ()4.68p . 

The eigenstates of Hq in the subspace N3 = n are labeled by subsets of size n of A<j = 
{1, 2, , . . . , s}; if K = {ki, k2, ■ ■ ■ , k^} is such a subset (where we will always assume 1 < A;i < 
k2 < ■ ■ ■ < kn < s), an eigenstate of Hq belonging to the eigenvalue 

n 

i=i 

is given by 

\Ek)= J2 V{K,M)\M). (5.2) 

MCAs;|M|=n. 
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For M = {xi , X2, • • • , Xn}, with 1 < xi < X2 < . . . < Xn ^ s, we have indicated above by | M) 
the simultaneous eigenstate of r3(l),r3(2) . . . ,T3(s) in which only the spins in M are "up", 
and we have set: 

V{K,M) = det{\\vkSx,)l,j=i,...,n) (5-3) 

where the functions have been defined in (|3.7|) . 
We set 

Qi\ {X1,X2, ■ ■ ■,Xn} ) = Xi\ Xi,X2, . . . ,X„ ). (5.4) 

It is easy to study, by the techniques of section [dl the asymptotic (as s — > +00 and t +00) 
joint distribution of the observables Qi, and therefore to give quantitative estimates of the 
correlation between the speeds of different particles and its dependence on the initial con- 
dition. To quote just one example, in the subspace N-^ = 2 and in the state | {1,2} ) the 
velocities {Vi,V2) of the two "up" spins (the limits in law of Qi/t and Q2/t, respectively) 
have joint probability density 

fVi,V2{vi,V2) = IiO,v2){n)I{0,l){V2) ^ r- ,,. . (5.5) 

It is immediate from (I5.5p to compute the conditional expectation -E'(T4|V2) of the velocity of 
the leftmost particle given the one of the rightmost particle; it turns out to be: 

E{Vi\V2) = ^ + 0{Vi). (5.6) 

In this section we advance the following idea: if the issue of the computation is the application, 
for a given number g of times, of a given primitive G to the register, initialize the cursor in 
the N3 = g subspace, in the state, say, \ {1,2, . . . , g} ); let then the system evolve according 
to the Hamiltonian: 

A '"^ 

H = - -'^Ucc(S)T+{x + 1)t_(x) + r+(x)T_(x + 1) (5.7) 

where 

= C'^o.-, for a fixed xq > 5, G° = I^. (5.8) 

An implementation of this approach is shown by the probability-entropy-Bloch diagram of 
figure [STTl Simple expressions for the quantities shown in figure [STT] can be obtained by the 
explicit form of the eigenvectors of the Hamiltonian described by (15. 7p and ()5.8p in every 
eigenspace of N3. For instance in the subspace = 3 a complete set of eigenstates is given, 
for C, = ±1 and l<j<h<k<s, by: 

lC;^{iAfe}) = ^ ^({i,/i, A;},{xi,X2,X3}) • 

l<Xl<X2<X-i<S 

. ^^(Xi-X0)+^(X2-X0)+^(X3-X0)| ^3 = ^)^1 {a;i,X2,X3} ) (5.9) 

where i} is the unit step function defined by: 
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1 r 



S(Pr(t)) 




(c) 

Figure 5.1: ^ — A, g — 3, = 6, s = 20, < t < 4s; Uxq = G = exp(— zq;(72/2) with a and 9 given 
by i4.44\ l o,d 14.4-3^ , Ux = Ir for x ^ xq; initial condition \ R{1) ) (X" | {1, 2, 3} ) with \ i?(l) ) given by 

!44W - 



5.2 Non commuting computational primitives 

The spectral structure ()5.9p is peculiar of the extremely simple situation ()5.8p (just one active 
link) considered there. As soon as we have more than one active link, say the primitive A 
acting on link (a, a + 1) and the primitive B acting on link {b,b + 1), with b > a + 1, a 
new phenomenon (that for simplicity we discuss in the = 2 case) takes place: the energy 
eigenstates have not anymore the form of a linear combinations of tensors products of the 
form M(xi,X2) \ CJ3 = C ) "X* | {^i, X2} ), with M{xi, X2) a monomial in A and B; related to this, 
the coordinates xi, X2 lose, strictly speaking, the meaning of relational time 47|: given that 
at a given value of t, Qi = xi and Q2 = X2 we can only claim that the state of the register 
has been acted upon by a polynomial in A and B. 

This phenomenon is easily understood in terms of the Dyson expansion of the propagator: 
the probability amplitude for the two excitations being in xi, X2 (both larger than b), given 
that at time they were in yi, 1/2 (both < a), receives contributions not only from Feynman 
paths along which the rightmost excitation goes past a and b and then the leftmost excitation 
goes past a and b (along such a computational path the state of the register is modified by 
BAB A), but also, among others, from paths along which both excitations go past a before 
both going past b (along such a computational path the state of the register is modified by 
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BBAA). 

By means of numerical simulation, it is possible to have some insight into the effect of the 
interaction by the register when it is acted upon by non-commuting operators. 
Let us consider a machine with one register spin in the state | 173 = +1 ) and a cursor chain 
of s sites initialized in the state | {1, 2} ), that is 

I ^^(0)) = I Ml) = 1^73 = +1)1 {1,2}). (5.11) 

Consider the following family of Hamiltonians 

s 

Ha,b = ^B^'''-A^''--T+{x + l)T^{x) + h.c., (5.12) 

x=l 

where 6j^x is the Kronecker 5 function. 

Each Hamiltonian of the family is characterized by having only two active links a — > a + 1 and 
5 — > & + 1, during which the operator A and B are respectively administered to the register. 
Figures 15.21 and 15.31 are suggestive of the following behavior: the motion of the excitations 
along the chain is unaffected by the interaction with the register as long as the active links are 
next to each other (b = a + 1). This reflects the fact that the energy spectrum of Hamiltonians 
of type Ha^a+i is the same as the one of the free Hamiltonian. In fact, it can be easily shown 
that the states 

lC;^b>}> = V({i,/i},{a;i,X2})- (5.13) 

l<a;i<X2<s 

where V({j, /i}) are defined as in ()5.3p . are eigenstates of Ha,a+i corresponding to the eigen- 
value Ey^iij defined in (|5.ip and 'd is defined as in (|5.10p . 

As soon as b > a + 1, namely the active links are separated by inactive links (that is links 
during which the identity is applied) localization appears. This suggests a change in the 
spectrum of Hamiltonians Ha^a+ki k > 1. Whereas in the first case (Ha^a+i) we can give 
analytic evidence of the equality of the "free particles" spectrum and the interacting one, in 
this second case we have only numerical evidence of what we stated. 
In the following examples, we choose the operators A = ai and B = (T3. We observe that 

aiasaiasl + 1 ) = e''^ aiaiasasl + 1 ) (5.14) 

Figures 15.21 15.31 and 15.41 show the time evolution, in the A''3 = 2 case, of the probability of 
reaching a determined configuration of the machine when b = a + 1 (solid black line) and 
b = a+2 (dashed red line) for different positions of the active links. For the first-neighbor case 
to all the computational paths ending in | g ) there corresponds the transformation BAB A 
on the register; in the second case, there are different computational paths ending in | g ) to 
which there correspond different transformations. When b = a + 2 the classes of paths for 
which the action on the register is BBAA = 1 and BAB A = — 1 have the same probability 
amplitude and gather the most of the probability mass. The destructive interference between 
them is then due to the change of phase (kickback effect, see (I5.14p ) of the amplitude brought 
about by the different transformations on the register. The consequent dynamical creation 
of boundary conditions causes a "confinement" of the traveling excitations on the left of the 
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Figure 5.2: s — 20, A — ai, B — the solid line represents the probability of observing 1 < Qi < 
Q2 < 5 when a — 2 and b — 3; the dashed line the same probability when a = 2 and 6 = 4. 




Figure 5.3: Same parameters and operators as in figure [KM We represent the probability Pt{Qi = 
9ij Q2 = 92) for t — A; (a) corresponds to the choice a = 2 and b — 3; (b) represents the same quantity 
for a = 2 and 6 = 4. 
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Pi(Ql=10, Q2 = 11, cr3=+1) P,(1<Qi<Q2< 11) 




(a) (b) 

Figure 5.^: s = 20, A = ai, B = a^; solid (black) line: a ~ 9, b ^ 10; dashed (red) line: a = 
9, 6 =11. (a) The destructive interference between the BBAA and BAB A causes a "reflection" 
on the dynamical boundary, (b) The "confinement" effect due to the destructive interference between 
different computational paths. 



interference location, here coinciding with the end of the last active link. 

Waiting for an algorithm that might benefit from the above possibility of simultaneously 
exploring different computational paths (concurrency ?), we explore, in figure [53} the idea 
(or classical prejudice?) that this nuisance can be in part avoided by using suitable initial 
conditions. The idea, suggested by ()4.70p . is of course to prepare on Aq = {1,2, ... ,o} an 
initial A^s = g state such that the excitations travel as spatially well localized wave packets 
of so different speeds that it is at any time unlikely that they simultaneously hit the region 
(a + 1,5). 

A remark about our insistence, throughout the previous chapter, in gathering experience 




Figure 5.5: /z = 4, s = 20, < i < 4s; (a) solid line: Ua = A, Ub = B, a = 6, 6 = 8, same initial 
condition as in fiaure \5.1[ (b) solid line: the = 3 state has been prepared by setting the initial chain 
{1, . . . , 6} in the ground state of the Hamiltonian X]x=i ''"+(^ + ^)^~ i^) + ^-^^ ^''^^ comparison purpose 
fiaure f5!71 a is reproduced in both frames as a dashed line. 



about the behavior of the evolution of a state of an initial subchain A^ = {1, 2, . . . , e} is now 
in order. 
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We observe that an initial state (not necessarily in the A^3 = 1 subspace) in of the form 



where Im is the indicator function of the set M and x = 1, . . . ,e, can be prepared as a 
post-kickback state (with respect to an ancilla qubit) after the reversible evaluation of a 
function / : {—1,1}'^ {—1,1}. We conjecture that subsequent evolution of | m ) under 
the Hamiltonian (j4.68p on = {1,2,... ,s}, with s » e, might help in setting tests of 
hypotheses about the Fourier coefficients 



of the function / via time-of- flight techniques. There is at least one non trivial case in which 
the above conjecture works: having prepared all spins in {e + l,...,s} in the "up" state, 
the Deutsch-Josza alternative [i^ " / constant (jcgl = 1) vs. / balanced (c0 = 0)" becomes 
equivalent to the alternative "stationary vs. non stationary" under the Hamiltonian ()4.68p . 
about the state of the overall system. 




(5.15) 




(5.16) 
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Chapter 6 

Conclusions and outlook 



The CCNOT gate can be implemented on Feynman's quantum computer by means of con- 
ditional jumps: the transition amplitudes at a vertex of the clocking graph depend on the 
state of the register. The Feynman machine, thus, not only anticipates the continuous time 
quantum walk paradigm but also extends it to what we call interacting quantum walks. Since 
the CCNOT alone constitutes a complete reversible logical basis, the computer model is uni- 
versal with respect to the class of function reversibly computable by a Turing machine. 

In dealing with specific algorithms, we abandoned the black-boxes, or oracular, setting and 
adopted a top-down approach. For instance in the case study of Grover's algorithm, we wrote 
down explicitly the Hamiltonian generating, at certain times, the evolution of the input state 
into the desired output state. In this specific context, we defined a scheme for the iteration 
of quantum subroutines and analyzed the space cost of it, which turns out to be logarithmic 
in the number of iteration. Moreover, we showed that the CCNOT can be implemented by 
using only three body interactions; this represents a significant reduction of the complexity of 
the physical realization of the computing device. 

We have shown that, by means of a generalized Peres basis, it is possible to express the evo- 
lution of the system as if the clocking excitation were traveling without interaction along a 
linear spin chain. As far as the dynamics of the interacting quantum walk is concerned, then, 
we can restrict our attention to the motion of the clocking excitation on a linear chain. 

The pure XY Hamiltonian Hq given in (14.680 describes, in the Luther-Liischer-Susskind for- 



malism [49|, |50|, [51[, a massless Dirac quantum field on a 1-dimensional lattice. The full 
Hamiltonian (11.20 is suggestive of the minimal coupling of this Fermi field, implementing 
the clock, with additional quantum fields implementing the register. We have been trying to 



contribute to the line of research, that seems to be emerging these days j52l. l53l. l54i|. devoted 
to making this connection between quantum computing and relativistic quantum field theory 
explicit. It is an easy guess that this quantum field theoretical intuition was well present in 
the original work [ij]. Particularly penetrating is, in this respect, Peres' remark that in Feyn- 
man's model calculations run forward and backward in time just as particles and antiparticles 



in Feynman's classical work on relativistic quantum field theory ( 18|], p. 3269). As a further 
remark, we observe that the three body interactions needed by Feynman's model are hard to 
conceive out of a field theoretical context. 

It is because of this field theoretical perspective that we have tried to avoid any "engineering" 
jHB] (space dependence) of the coupling constant A in (14.680 . well aware of the fact that, in the 
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Dirac XY correspondence, A is related to the spacing adopted in the lattice approxima- 
tion. In such a context it would be very hard to understand (without a projection mechanism 



551], which seems to have an exponential cost) the implementation of a space dependence 



such as 



X{x) = const. \J x{x — s) (6-1) 



that leads in and 55| to the existence of sharply distinguished instants in which the 
position of the cursor is certain. Nor would it be easy to understand (j6.ip in a solid state 
implementation [s^, where A is related to the effective mass of the clocking excitation. 

We have focused our attention on the clocking field t{x), singled out as the one which, un- 
der suitable boundary conditions and for initial conditions localized close to the boundary, 
exhibits particle-like excitations performing, for long enough intervals of t, a quantum walk 
in a distinguished direction. 

Spatial homogeneity of the chain leads to the existence of the limit in law V = limt^+oo Q{t)/t 
for the position of such an excitation on a semi-infinite {s — > +oo) box. In the A^3 = 1 sub- 



space, because of Peres' conservation law [18|, the observable Q acquires the meaning of rela- 
tional time (given the observed value of Q, the state of the register is known with certainty) 
and, therefore, the random variable V acquires the meaning of number of computational steps 
per unit t. The fact that the variance of V is strictly positive has the effect that in terms of 
the parameter time t (as opposed to relational time Q) the evolution of the register appears 
to be dissipative: we have, for a simple model, written the corresponding Lindblad evolution 
and studied the ensuing buildup of entropy. 

On our simple instance of quantum search we have shown that, in the "low level", physical 
approach that we pursue (in which time runs, for the register, because it is coupled with an 
additional quantum field) the buildup of entropy imposes an upper bound on the probability 
of finding the target state which is more severe than the one predicted by the "high level", 
algorithmic approach (in which the successive primitives are applied by an external macro- 
scopic agent). 

In the attempt of decreasing the deficit in the probability of success in a quantum search, due 
to the decohering effect of the coupling with the clocking field, we have provided examples 
of the benefit of spreading the initial wave function of the cursor on an initial launch pad 
instead of, as a classical prejudice would suggest, having it strictly localized at one site. 

The master equation of the state of the register satisfies the Lindblad equation; therefore, if 
on one side the clock makes the register evolve, on the other it acts as a decohering environ- 
ment. 

We point out, furthermore, that the analysis of the probability cost of an algorithm will be 
significant only when a optimal preparation-measurement scheme will be found. The n-pulse 
trap and the launch-pad are proposals in this direction. We observe that, with the same 
amount of space resources, a longer sojourn time can be achieved by means of an iteration of 
quantum subroutines circuit in which to trap the excitation. 

The study of the energy distribution of the system before and after the (non-optimal) mea- 
surement allows for the analysis of the resetting cost of the machine which, as long as we 
consider perfect chains, is the only energetic cost of the computation. 

The role of initial condition must be studied also when the initial state of the machine presents 
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entanglement between the register and the clock; for example let's take the initial condition 

s d 
x=l j=l 

with {| rj )}j=i is a selected orthonormal basis of Hregister ■ 

With the same techniques used in section 14.11 it is possible to show that 

cPv(io)i^) = lim Urn (eo|e^*''expfiz^')e-**^Uo) = 

i2i;l^(j>rcsin(t;))p + vr - arcsin(u))| 



1=1 



\/l - ■v^ 



-dv (6.3) 



where 



^(^') = E( 1^1^ • • • UU\ n )i{i,x) (6.4) 

X = l 1 = 1 

is the sine transform of the initial state. The interest of the most general form of the initial 
vector state relies on the possibility of exploiting the entanglement between the register and 
the cursor in order to either increase the average speed of computation or reduce the spreading 
of the wave packet. The analytic description of the interaction between the two subsystems 
models also some kind of decoherence appearing in spin chains. 

We have, moreover, abandoned the classical prejudice of having a single clocking excitation, 
providing a preliminary analysis of the idea of starting the cursor in an initial state with 
iVs > 1. We have shown, in this context, an efficient way of iterating the application of a 
single primitive to the register. Numerical examples of section 15.21 show that, however, as 
soon as the unitary operators acting on the register do not commute, the positions of the 
excitations have no longer the meaning of relational time of the system: the state of the reg- 
ister depends also on the path made by the excitations along the chain. Indeed, we showed 
numerical examples in which interference between different paths leads to localization a la 
Anderson of the excitations. In that case, however, it is the presence of impurities in the 



lattice, or chain, which leads to localized states. We refer to 57| for a recent discussion of 
Anderson localization on spin chains. 

The case > 1 deserves, we think, further research, both from the algorithmic and the 
physical point of view. 

From the algorithmic point of view we plan to examine other instances (beyond the one cur- 
sorily examined at the end of section 15. 2p in which time-of-flight spectroscopy (based on the 
Fourier transform vs. speed relationship recalled in section 14. ip of the post-kickback state 
can answer Yes/No questions about the algorithm. 

Prom the physical point of view, the "obvious" choice of the "all down" reference state made 
throughout this work is far from being optimal from the point of view of studying the thermo- 
dynamic cost of resetting the register. The best reference state for the study of this ultimate 
cost of reversible computation would of course be the ground state and, for Hamiltonians 
of the form (j4.68p . with s even, it is an = s/2 state. This will require, we think, the 
formulation of an appropriate Bethe Anzatz for the Hamiltonian ()1.2p . 
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It is this last point which leads to some intrinsic difficulty in the analytic description of the 
system. In fact, if scattering theory includes an analytic treatment of the perturbed Hamil- 
tonian for the XY chain, it seems quite hard to do the same if the coefficients of the matrix 
belong to a non-Abelian algebra, which is our general case. Nevertheless, we have shown that 
the interaction between the register and the multi-handed clock shows new effects such as the 
confinement of the excitation resembling Anderson localization. 

In section [5^ we have related the probability deficit of figure [53l a) to the simultaneous pres- 
ence of two of the excitations in the region between the active links a — > a -|- 1 and b ^ b + 1, 
which corresponds to the application in the "wrong" order of the oracle and estimation step 
of Grover's algorithm. If we initialize the initial condition of the cursor in the ground state 
of the launch pad {1, 2, . . . , 6} we partially fill the probability gap due to the sojourn of the 
excitations in the critical region. Similarly, the confinement effect shown in figures 15.21 and 
15.41 is related to the interference pattern between states of the form 

BABA\as = +l)\{qi,q2}) = e'^ks = +1 )| M, ^2} ) (6.5) 

and states of the form 

BBAA\a3 = +l)\{qi,q2}) = ks = +1 )| {gi, 92} ), (6.6) 

with b + 1 < qi < q2 < s. 

The effect is similar to scattering on an impurity: the relative phase e"^ between (j6.5p and 
(16. 6p leads, by a familiar kickback mechanism in which the register plays the role of an ancilla 
qubit, to destructive interference in the distribution of the cursor. 

The transmission amplitude is proportional to {a — f3), a being the amplitude associated to 
states of the form ()6.5p and (3 the amplitude of states of the form (j6.6p . 

The trasmission probability deficit is compatible with the empirical distribution of the random 
variable sojourn time in the critical region {(xi, X2) : a < xi < X2 < b}, with a = 9 and b = 11, 
shown in figure [6?T1 
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Figure 6.1: s ~ 20; H = Hq as in {4^.6^ ; — 2; a ^ 9, h ~ 11; initial condition \ {1, 2} \ 
1000 paths. The empirical distribution function of the sojourn time in the critical region. 



sample of 



We define (see ^]) the sojourn time of a quantum walk in a given region as the sojourn time 
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of a canonically associated controlled random walk. The idea that "stochastic control theory 
can provide a very simple model simulating quantum mechanical behavior" is borrowed from 
59( 1 ■ The actual details of the numerical simulation in our discrete settings implements the 
prescription of [60]. 

What figure ET] says is that 80% of the trajectories of our sample, of size 1000, behave as the 
one shown in figure [621 (a), never hitting the vertex (11, 10); only 20% have at some time both 
excitations between a + 1 and b. As a sketch of the body of ideas relating quantum walks 
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Figure 6.2: Same parameters as in fiaure \KT[ (a) A trajectory not passing through the critical region, 
(h) A trajectory traversing the critical region. The critical region, that is the configuration (11, 10), is 
highlighted by a (green) point. 



to random walks [61[ we discuss the example referring to the case = 1 represented in 
figure [631 The whole point behind the algorithm leading to figure [6T3l is that the probability 
density 



4x^ 

p(t,x) = -^J^{t) 
corresponding to the amplitude (j3.13p satisfies the equation 



(6.7) 



d , , 



-sign 



{JS)Jx+l{t))yJp{t, x)p{t, X + l) + 



+ sign{J^{t)Jx-i{t))yJ p{t,x - l)p{t,x). 



(6.8) 



Equation (j6.8p can, in turn, be read as the continuity equation for a birth and death process 
on Aoo = {1, 2, . . .}. The sample paths (t, q{t)) of this process are shown in figure [Ol 

The stochastic simulation of quantum phenomena discussed in [s^, [H^, [60I raises extremely 
interesting questions about the foundations of quantum mechanics itself [64]. 

We wish to stress as a final remark that stochastic mechanics poses also interesting questions 
in the theory of quantum computation: is the exponential speedup offered by quantum walks 
in crossing a graph or decision tree 19|, [2l| attainable by a classical stochastic algorithm? 
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